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1  Introduction 


Objectives  and  Scope 

Current  safety  analysis  procedures  and  computer  programs  such  as  CSLIDE 
(Pace  and  Noddin  1987)  for  concrete  retaining  walls  and  gravity  structures  com¬ 
pute  traditional  safety  factors  that  are  not  based  on  reliability  analysis.  CSLIDE 
defines  the  safety  factor  as  the  ratio  of  sliding  resistance  to  sliding  force.  Safety 
factors  are  not  accurate  measures  of  stability  reliability  because  they  do  not 
account  for  the  various  uncertainties  in  underlying  parameters  or  variables  of 
stability  problems.  Also,  the  factors  of  safety  do  not  convey  the  nonlinear  nature 
of  relationships  between  the  margin  of  safety  that  they  measure  and  the  unsatis¬ 
factory  performance  likelihood  that  can  be  used  as  a  basis  for  measuring  stability 
reliability. 

The  objective  of  this  study  was  to  develop  reliability  assessment  methods  for 
the  stability  of  gravity  concrete  structures.  This  report  describes  the  probability- 
based  reliability  assessment  methodology  for  concrete  retaining  walls  and  gravity 
structures  that  was  developed.  The  methodology  is  based  on  the  U.S.  Army  Corps 
of  Engineers  Computer-Aided  Structural  Engineering  (CASE)  program  for  Sliding 
Stability  of  Concrete  Structures  (CSLIDE,  Pace  and  Noddin  1987).  A  user 
interface  for  CSLIDE  and  the  reliability  program  based  on  CSLIDE  (called 
RCSLIDE)  was  also  developed  using  Microsoft  Visual  Basic.  The  sigiiificance  of 
the  software  development  procedure  described  in  this  study  is  that  is  establishes 
prototype  reliability  software  that  is  modular  and  based  on  an  existing  CASE 
program.  Other  CASE  programs  can  be  modified  and  utilized  for  reliability 
purposes  in  a  similar  fashion. 

The  development  of  the  methodology  required  the  definition  of  a  performance 
function  for  retaining  walls  and  gravity  structures,  development  of  a  library  of 
probabihty  functions  (Ayyub  and  Chao  1994),  development  of  a  structural  reli¬ 
ability  assessment  module,  development  of  user  interfaces,  and  selection  and 
performance  of  test  cases. 
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Report  Structure 

Chapter  2  describes  numerical  algorithms  and  basic  theory  for  the  following 
two  probability-based  reliability  assessment  methods:  (a)  advanced  second 
moment  (ASM)  method,  and  (b)  Monte  Carlo  simulation  (MCS)  using  direct  (DR) 
or  importance  sampling  (IS).  Chapter  3  describes  the  stability  analysis  of  concrete 
gravity  structures  and  retaining  walls  and  the  CASE  computer  program  CSLIDE. 
In  Chapter  4,  numerical  algorithms  for  nonclosed  form  performance  function  used 
in  the  reliability  analysis  methods  and  the  computer  software  RCSLIDE  (relia¬ 
bility  and  stability  assessment  of  concrete  gravity  structures)  are  developed.  In 
Chapter  5,  two  case  studies  with  11  random  variables  are  presented  to  illustrate 
the  use  of  the  developed  reliability  assessment  software  RCSLIDE.  The  case 
studies  deal  with  noncorrelated  and  correlated  random  variables.  Conclusions  and 
recommended  future  work  are  provided  in  Chapters  6  and  7,  respectively. 


Companion  Manual 

The  CSLIDE  manual  (Pace  and  Noddin  1987)  referenced  in  this  report  is  in 
the  process  of  being  updated  to  reflect  the  Windows  version  of  the  program  and 
the  incorporation  of  reliability  assessment.  This  new  report,  “Reliability  and 
Stability  Assessment  of  Concrete  Gravity  Structures:  User’s  Guide,”  will  also  be 
a  U.S.  Army  Engineer  Waterways  Experiment  Station  technical  report  and  a 
companion  manual  to  this  theoretical  manual. 
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2  Structural  Reliability 
Assessment 


The  reliability  of  an  engineering  system  can  be  defined  as  its  ability  to  fulfill 
its  design  purpose  for  some  time  period.  The  theory  of  probability  provides  the 
fundamental  basis  to  measure  this  ability.  The  reliability  of  a  structure  can  be 
viewed  as  the  probability  of  its  satisfactory  performance  according  to  some  per¬ 
formance  functions  under  specific  service  and  extreme  conditions  within  a  stated 
time  period.  In  estimating  this  probability,  system  uncertainties  are  modeled  using 
random  variables  with  mean  values,  variances,  and  probability  distribution  func¬ 
tions.  Many  methods  have  been  proposed  for  structural  reliability  assessment 
purposes,  such  as  first-order  second  moment  (FOSM)  method,  ASM  method,  and 
computer-based  MCS  (e.g.,  Ang  and  Tang  1990,  Ayyub  and  Haidar  1984,  White 
and  Ayyub  1985,  Ayyub  and  McCuen  1997).  In  this  chapter,  two  probabilistic 
methods  for  reliability  assessment  are  described:  (a)  the  ASM  method,  and  (b)  the 
MCS  method  using  DR  and  IS. 


Advanced  Second  Moment  (ASM)  Method 

The  reliability  of  a  stmcture  can  be  determined  based  on  a  performance  func¬ 
tion  that  can  be  expressed  in  terms  of  basic  random  variables  Xi’s  for  relevant 
loads  and  stmctural  strength.  Mathematically,  the  performance  function  Z  can  be 
described  as 

Z  =  Z ,  X2 , . . . ,  )  =  Structural  strength  -  Load  effect  (1) 

where  Z  is  called  the  performance  function  of  interest.  The  unsatisfactory  per¬ 
formance  surface  (or  the  limit  state)  of  interest  can  be  defined  as  Z  =  0.  Accord¬ 
ingly,  when  Z  <  0,  the  structure  is  in  the  unsatisfactory  performance  state,  and 
when  Z  >  0,  it  is  in  the  safe  state.  If  the  joint  probability  density  function  for  the 
basic  random  variables  Z,-  ’s  is  X2,...,x  (^i’  unsatisfactory 

performance  probability  7^  of  a  structure  can  be  given  by  the  integral 
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(2) 


Pu  ~f"'f  fXi,X2,...,X„  (Xi,X2,---,Xfi)dXi  dX2  ...dx„ 

where  the  integration  is  performed  over  the  region  in  which  Z  <  0.  In  general, 
the  joint  probability  density  function  is  unknown,  and  the  integral  is  a  formidable 
task.  For  practical  purposes,  alternate  methods  of  evaluating  are  necessary. 


Reliability  index 

Instead  of  using  direct  integration  as  given  by  Equation  2,  the  performance 
function  Z  in  Equation  1  can  be  expanded  using  a  Taylor  series  about  the  mean 
value  of  Z’s  and  then  tmncated  at  the  linear  terms.  Therefore,  the  first-order 
approximate  mean  and  variance  of  Z  can  be  shown,  respectively,  as 


(3) 


and 


a 


2 

Z 


n  n 


(dZ\ 

dZ  \ 

[dX,) 

dXi 
^  / 

Cov(Xi,Xj) 


(4) 


where 


=  mean  of  Z 

[X  =  mean  of  a  random  variable 
a|  =  variance  of  Z 

Cov(XijX j)  =  the  covariance  of  X^  and  X j 

The  partial  derivatives  of  dZjdXi  are  evaluated  at  the  mean  values  of  the  basic 
random  variables.  For  statistically  independent  random  variables,  the  variance 
expression  can  be  simplified  as 


(5) 


A  measure  of  reliability  can  be  estimated  by  introducing  the  reliability  index  p 
that  is  based  on  the  mean  and  standard  deviation  of  Z  as 
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p 


(6) 


Oz 

If  Z  is  assumed  to  be  normally  distributed,  then  it  can  be  shown  that  the  unsatis¬ 
factory  performance  probability  is 

P„=l-a»(P)  (7) 


where  <I>  is  the  cumulative  distribution  function  of  standard  normal  variate. 

The  aforementioned  procedure  of  Equations  3  to  7  produces  accurate  results 
when  the  performance  function  Z  is  normally  distributed  and  linear. 


Nonlinear  performance  functions 

For  nonlinear  performance  functions,  the  Taylor  series  expansion  of  Z  is 
linearized  at  some  point  on  the  unsatisfactory  performance  surface  called  design 
point,  checking  point,  or  the  most  likely  unsatisfactory  performance  point  rather 
than  at  the  mean.  Assuming  the  original  basic  variables  Xi  ’s  are  uncorrelated, 
the  following  transformation  can  be  used: 


~  M' jf,. 


(8) 


If  Xf ’s  are  correlated,  they  need  to  be  transformed  to  uncorrelated  random  vari¬ 
ables,  as  described  by  Thoft-Christensen  and  Baker  (1982),  Ayyub  and  McCuen 
(1997),  or  Ang  and  Tang  (1990)  and  described  in  the  later  section  entitled 
“Correlated  random  variables.”  The  reliability  index  P  is  defined  as  the  shortest 
distance  to  the  unsatisfactory  performance  surface  from  the  origin  in  the  reduced 
y-coordinate  system.  The  point  on  the  unsatisfactory  performance  surface  that 
corresponds  to  the  shortest  distance  is  the  most  likely  unsatisfactory  performance 
point.  Using  the  original  Z-coordinate  system,  the  reliability  index  P  and  design 

point  (Zj* ,  Zj , . . . ,  Z  * )  can  be  determined  by  solving  the  following  system  of 
nonlinear  equations  iteratively  for  P : 


( bz 

dX, 


il/2 


vf— 


'at, 


(9) 
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(10) 


z(x;,z2,...,z;)=o  (11) 

where  a,'  is  a  directional  cosine  and  the  partial  directives  are  evaluated  at  the 
design  point,  indicated  by  an  asterisk.  Then,  Equation  7  can  be  used  to  evaluate 
.  However,  this  formulation  is  hmited  to  normally  distributed  random  vari¬ 
ables.  In  reliability  assessment,  the  directional  cosines  can  be  viewed  as  measures 
of  the  importance  of  the  corresponding  random  variables  in  determining  the 
reliability  index  p.  Also,  partial  safety  factors  y  that  are  used  in  load  and 
resistance  factor  design  (LRFD)  can  be  computed  as  follows: 


In  general,  partial  safety  factors  take  on  values  larger  than  one  for  load  variables 
(called  in  this  case  load  amplification  factors),  and  values  less  than  one  for 
strength  variables  (called  in  this  case  strength  reduction  factors). 


Equivalent  normal  distributions 

If  a  random  variable  X  is  not  normally  distributed,  then  it  needs  to  be  trans¬ 
formed  to  an  equivalent  normally  distributed  random  variable,  indicated  by  the 

superscript  N.  The  parameters  of  the  equivalent  normal  distribution,  and 

of ,  can  be  estimated  by  imposing  two  conditions  (Rackwitz  and  Fiessler  1976, 
1978).  The  cumulative  distribution  functions  (CDF)  and  probability  density  func¬ 
tions  of  a  nonnormal  random  variable  and  its  equivalent  normal  variable  should  be 
equal  at  the  design  point  on  the  unsatisfactory  performance  surface.  The  first 
condition  can  be  expressed  as 


O 


^Xi 


N  \ 


The  second  condition  is 


(13) 


^Xi 


N  \ 


(14) 
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where 


Fj  =  nonnormal  cumulative  distribution  function 

(j)  =  probability  density  function  of  standard  normal  variate 

f.  =  noimormal  probability  density  function 


The  standard  deviation  and  mean  of  equivalent  normal  distributions  can  be  shown, 
respectively,  to  be 


<  = 


and 


(15) 


(16) 


Having  determined  a  ^  and  p  ^  for  each  random  variable,  p  can  be  solved  using 

A  i  A I 

the  same  procedure  of  Equations  9-11. 

The  ASM  method  is  capable  of  dealing  with  nonlinear  performance  functions 
and  noimormal  probability  distributions.  However,  the  accuracy  of  the  solution 
and  the  convergence  of  the  procedure  depend  on  the  nonlinearity  of  the  perform¬ 
ance  function  in  the  vicinity  of  design  point  and  the  origin.  If  there  are  several 
local  minimum  distances  to  the  origin,  the  solution  process  may  not  converge  onto 
the  global  minimum.  The  probability  of  unsatisfactory  performance  is  calculated 
from  the  reliability  index  P  using  Equation  7,  which  is  based  on  normally  distrib¬ 
uted  performance  functions.  Therefore,  the  resulting  unsatisfactory  performance 
probability  based  on  the  ASM  is  approximate  except  for  linear  performance 
functions  because  it  does  not  account  for  any  nonlinearity  in  the  performance 
functions. 


Correlated  random  variables 

Reliability  analysis  of  gravity  structures  needs  to  be  based  on  correlated  soil 
properties  such  as  angle  of  internal  friction  and  cohesion  for  soil  layers.  In  this 
section,  this  correlation  is  assumed  to  occur  between  pairs  of  random  variables  for 
each  layer.  Also,  correlated  random  variables  are  assumed  to  be  normally  distrib¬ 
uted  since  nonnormal  and  correlated  random  variables  require  additional  informa¬ 
tion  such  as  their  joint  probabihty  density  function  or  conditional  distributions  for 
their  unique  and  full  definition  (Ang  and  Tang  1990).  Such  information  is  com¬ 
monly  not  available  and  difficult  to  obtain.  A  correlated  (and  normal)  pair  of 
random  variables  Xi  andX2  with  a  correlation  coefficient  p  can  be  transformed 
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into  noncorrelated  pair  Yi  and  Y2  by  solving  for  two  eigenvalues  k  and  the 
corresponding  eigenvectors  as  follows: 


2t 


(X 


\ 


1  ~  M-Xi 


+ 


X2~  1^X2 
0X2 


\ 


/ 


(17) 


1  /  Xi  - 

2t  ^  Xi 

\ 


X2  ~  '' 


a, 


(18) 


where  t  =  VCL5  .  The  resulting  Y  variables  are  noncorrelated  with  respective  vari¬ 
ances  that  are  equal  to  the  eigenvalues  k  as  follows: 


Ojr  —  X,j  —  1  +  p 


(19) 


ay^=X2=l-p  (20) 

For  a  correlated  pair  of  random  variables,  Equations  9  and  10  need  to  be  revised 
respectively  to  the  following: 


(23) 

X2  =  1^X2  ~  ) 

(24) 

where  the  partial  derivatives  are  evaluated  at  the  design  point. 
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Numerical  algorithms 


The  ASM  method  can  be  used  to  assess  the  rehability  of  a  structure  according 
to  a  nonhnear  performance  function  that  may  include  nonnormal  random  vari¬ 
ables.  Also,  the  performance  function  can  be  in  a  closed  or  nonclosed  form 
expression.  The  implementation  of  this  method  requires  the  use  of  efficient  and 
accurate  numerical  algorithms  in  order  to  deal  with  the  nonclosed  forms  for 
performance  function.  The  ASM  algorithm  can  be  summarized  by  the  following 
steps  using  two  cases: 

Algorithm  2-1  (noncorrelated  random  variables).  Use  the  following  steps: 

(1)  Assign  the  mean  value  for  each  random  variable  as  a  starting  design  point 
value,  i.e., 

(2)  Compute  the  standard  deviation  and  mean  of  the  equivalent  normal  distri¬ 
bution  for  each  nonnormal  random  variable  using  Equations  13-16. 

(3)  Compute  the  partial  derivative  dZjdXi  of  the  performance  function  with 
respect  to  each  random  variable  evaluated  at  the  design  point  as  needed  by 
Equation  9. 

(4)  Compute  the  directional  cosine  a,  for  each  random  variable  as  given  in 
Equation  9  at  the  design  point. 

(5)  Compute  the  rehability  index  P  by  substituting  Equation  10  into 
Equation  11  and  satisfying  the  hmit  state  Z  =  0  in  Equation  11  using  a 
numerical  root-finding  method. 

(6)  Compute  a  new  estimate  of  the  design  point  by  substituting  the  resulting 
rehability  index  P  obtained  in  Step  5  into  Equation  10. 

(7)  Repeat  Steps  2  to  6  until  the  reliability  index  p  converges  within  an 
acceptable  tolerance  6. 

End  of  Algorithm  2-1. 

Algorithm  2-2  (correlated  random  variables).  Use  the  following  steps: 

(1)  Assign  the  mean  value  for  each  random  variable  as  a  starting  design  point 

value,  i.e.,  {Xl,X*2,...,xl)  =  , ^^2 > •  •  • . )  • 

(2)  Compute  the  standard  deviation  and  mean  of  the  equivalent  normal  distri¬ 
bution  for  each  nonnormal  random  variable  using  Equations  13-16. 

(3)  Compute  the  partial  derivahve  dZfdX^  of  the  performance  function  with 
respect  to  each  noncorrelated  random  variable  evaluated  at  the  design 
point  as  needed  by  Equation  9. 
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(4)  Compute  the  directional  cosine  a,  for  each  noncorrelated  random  variable 
as  given  in  Equation  9  at  the  design  point.  For  correlated  pairs  of  random 
variables,  Equations  21  and  22  should  be  used. 

(5)  Compute  the  rehability  index  P  by  substituting  Equations  10  (for  noncor¬ 
related  random  variables)  and  23  and  24  (for  correlated  random  variables) 
into  Equation  11  and  satisfying  the  limit  state  Z  =  0  in  Equation  11  using 
a  numerical  root-finding  method. 

(6)  Compute  a  new  estimate  of  the  design  point  by  substituting  the  resulting 
rehability  index  P  obtained  in  Step  5  into  Equations  10  (for  noncorrelated 
random  variables)  and  23  and  24  (for  correlated  random  variables). 

(7)  Repeat  Steps  2  to  6  until  the  reliability  index  p  converges  within  an 
acceptable  tolerance  6. 

End  of  Algorithm  2-2. 


Monte  Carlo  Simulation  (MCS)  Methods 

Monte  Carlo  simulation  techniques  are  basically  samphng  processes  that  are 
used  to  estimate  the  unsatisfactory  performance  probability  of  a  stmcture.  The 
basic  random  variables  in  Equation  1  are  randomly  generated  (Press  et  al.  1992; 
Ayyub  and  McCuen  1997;  Ayyub  and  Chao  1994)  and  substituted  into  Equa¬ 
tion  1.  Then  the  fraction  of  cases  that  resulted  in  unsatisfactory  performance  are 
determined  in  order  to  assess  unsatisfactory  performance  probability.  Two 
methods  are  described  in  this  section:  the  direct  MCS  and  the  IS. 


Direct  MCS  method 

The  direct  simulation  method  comprises  drawing  samples  of  the  basic  noncor¬ 
related  variables  according  to  their  corresponding  probabilistic  characteristics  and 
then  feeding  them  into  the  performance  function  Z  as  given  by  Equation  1. 
Assuming  Nu  to  be  the  number  of  simulation  cycles  for  which  Z  <  0  in  a  total  N 
simulation  cycles,  then  an  estimate  of  the  mean  unsatisfactory  performance 
probability  can  be  expressed  as 


P 

“  N 


(25) 


The  estimated  unsatisfactory  performance  probability  P„  should  approach  the 
true  value  for  the  population  when  N  approaches  infinity.  The  variance  of  the 
estimated  unsatisfactory  performance  probability  can  be  approximately  computed 
using  the  variance  expression  for  a  binomial  distribution  as 
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Fa/-(P„)  = 


i^-Pu)Pu 

N 


Therefore,  the  coefficient  of  variation  (COV)  of  the  estimated  unsatisfactory 
performance  probability  is 


(27) 

These  equations  show  that  direct  simulation  can  be  economically  prohibitive  in 
some  cases,  especially  for  small  unsatisfactory  performance  probabilities.  In  the 
subsequent  section,  the  IS  method  is  described  for  the  purpose  of  increasing  the 
efficiency  of  simulation. 


Importance  sampling 

The  probability  of  unsatisfactory  performance  of  a  structure  according  to  the 
performance  function  of  Equation  1  is  provided  by  the  integral  of  Equation  2.  In 
the  evaluation  of  this  integral  with  direct  simulation,  the  efficiency  of  the  simula¬ 
tion  process  depends  on  the  magnitude  of  the  probability  of  unsatisfactory  per¬ 
formance,  i.e.,  the  location  of  the  most  likely  unsatisfactory  performance  point  or 
design  point  as  defined  in  the  section  “Nonlinear  performance  functions.”  The 
larger  the  margin  of  safety  (Z)  and  the  smaller  its  variance,  the  larger  the  simula¬ 
tion  effort  needed  to  obtain  sufficient  simulation  runs  with  unsatisfactory  perform¬ 
ances;  in  other  words,  smaller  unsatisfactory  performance  probabilities  require 
larger  numbers  of  simulation  cycles.  This  deficiency  can  be  addressed  by  using 
importance  sampling.  In  this  method,  the  basic  random  variables  are  generated 
according  to  some  carefully  selected  probability  distributions  (called  importance 
density  function,  with  mean  values  that  are  closer  to  the  design  point  than 
their  original  (actual)  probability  distributions.  It  should  be  noted  that  the  design 
point  is  not  known  in  advance;  the  analyst  can  only  guess.  Therefore,  simulation 
runs  with  unsatisfactory  performances  are  obtained  more  fi-equently  and  the  simu¬ 
lation  efficiency  is  increased.  To  compensate  for  the  change  in  the  probability 
distributions,  the  results  of  the  simulation  cycles  should  be  corrected.  The  funda¬ 
mental  equation  for  this  method  is  given  by 

--y/, 

where 

/  =  performance  indicator  function  that  takes  values  of  either  0 
for  unsatisfactory  performance  or  1  for  survival 
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f^Xu,X2i, x„,)  =  the  original  joint  density  function  of  the  basic  random 

variables  evaluated  at  the  f  generated  values  of  the  basic 
random  variables 

h^xii,X2i, x„i)  =  the  selected  joint  density  function  of  the  basic  random 

variables  evaluated  at  the  f'  generated  values  of  the  basic 
random  variables 

For  noncorrelated  basic  random  variables,  the  joint  density  function ^A:ii,X26 
x„i)  can  be  replaced  by  the  product  of  the  density  functions  of  the  individual 
random  variables.  Similarly,  the  joint  h^Xii,X2b  Xni)  can  be  replaced  by  the 
product  of  the  corresponding  importance  density  functions.  In  Equation  28,  %  (x) 
is  called  the  sampling  (or  weighting)  density  function  or  the  importance  function. 
Efficiency  (and  tos  the  required  number  of  simulation  cycles)  depends  on  the 
choice  of  this  sampling  density  function.  The  coefficient  of  variation  of  the  esti¬ 
mated  unsatisfactory  performance  probability  can  be  based  on  the  variance  of  a 
sample  mean  as  follows: 


C07(P„)  = 

Additional  information  on  importance  sampling  can  be  found  in  SchueUer  and  Stix 
(1987). 

Correlated  random  variables 

In  this  section,  correlation  between  pairs  of  random  variables  is  treated  for 
simulation  purposes.  Similar  to  the  previous  section  in  this  chapter  on  correlated 
random  variables  for  ASM,  correlated  random  variables  are  assumed  to  be 
normally  distributed  since  nonnormal  and  correlated  random  variables  require 
additional  information  such  as  marginal  probability  distribution  for  their  unique 
and  full  definition.  Such  information  is  commonly  not  available  and  difficult  to 
obtain.  A  correlated  (and  normal)  pair  of  random  variables  Xi  and  X2  with  a 
correlation  coefficient  p  can  transformed  using  linear  regression  transformation  as 
follows: 


1  ^ 

1  V 

’  ^2i  ^  ^ni  ) 

X2  “  +  b-yX-^  +  8 


(30) 


where 

bo  =  the  intercept  of  a  regression  line  between  Xi  and  X2 
bi  =  slope  of  the  regression  line 
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s=  random  (standard)  error  with  a  mean  of  zero  and  a  standard  deviation  as 
given  in  Equation  33 

These  regression  model  parameters  can  be  determined  in  terms  of  the  probabilistic 
characteristics  ofXi  and-Y2  as  follows; 


(31) 

bo  =  1^X2  “  bi\ix^ 

(32) 

(33) 

The  simulation  procedure  for  a  correlated  pair  of  random  variables  (Xi  and  X2) 
can  then  be  summarized  as  follows: 

a.  Step  1:  Compute  the  intercept  of  a  regression  line  between  andX2  (bo), 
the  slope  of  the  regression  line  (bi),  and  the  standard  deviation  of  the 
random  (standard)  error  (e)  using  Equations  31-33. 

b.  Step  2:  Generate  a  random  (standard)  error  using  a  zero  mean  and  a 
standard  deviation  as  given  by  Equation  33. 

c.  Step  3:  Generate  a  random  value  forXi  using  its  probabilistic  charac¬ 
teristics  (mean  and  variance). 

d.  Step  4:  Compute  the  corresponding  value  of  X2  as  follows  (based  on 
Equation  30): 

X2=bQ+  bpc^  +  E  (34) 

where  bo  and  bi  are  computed  in  Step  1;  eis  a  generated  random 
(standard)  error  from  Step  2;  and  x\  is  a  generated  value  from  Step  3. 

e.  Step  5:  Use  the  resulting  random  (but  correlated)  values  of  xi  and  X2  in  the 
simulation-based  reliability  assessment  methods. 

This  procedure  is  applicable  for  both  the  DR  method  and  IS.  In  the  case  of  IS, 
correlated  random  variables  should  not  be  selected  for  defining  the  sampling  (or 
importance)  density  function  (h£). 
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3  Stability  Analysis  of 
Concrete  Gravity 
Structures 


The  reliability  assessment  of  the  sliding  stability  of  gravity  concrete  structures 
requires  the  definition  of  strength,  loads,  and  reliabihty  methods.  In  this  chapter, 
the  sliding  stability  strength  and  loads  needed  for  reliability  assessment  are  briefly 
described.  The  strength  and  loads  models  were  based  on  the  models  that  are  used 
in  the  sliding  stability  for  concrete  structures  (CSLIDE)  program  (Pace  and 
Noddin  1987). 

The  first  section  of  this  chapter  provides  theoretical  background  information 
on  computing  the  strength  and  load  for  the  sliding  stability  for  a  concrete  gravity 
structure.  The  only  failure  mode  that  is  handled  in  this  study  is  the  sliding 
stability  mode.  In  the  second  section,  the  CSLIDE  computer  program  is  briefly 
described  and  referenced. 


Theoretical  Basis 

Gravity  concrete  structures  serve  a  vital  function  in  many  Corps  facilities  such 
as  navigation  locks  and  dams.  A  gravity  structure  supports  earth  and  water  pres¬ 
sures  on  both  its  sides  and  bottom.  The  purpose  of  sliding  stability  analysis  is  to 
assess  the  ability  of  the  structure  in  fulfilling  its  function  without  an  excessive 
horizontal  deformation.  An  excessive  horizontal  deformation  can  be  viewed  as  a 
unsatisfactory  performance  of  the  structure.  In  order  to  make  such  an  assessment, 
the  sliding  strength  or  resistance  in  the  form  of  its  maximum  shear  strength  R  and 
the  applied  load  L  in  the  form  of  an  applied  shear  need  to  be  evaluated.  A  safety 
factor  SF  or  margin  of  safety  can  then  be  estimated  as 

SF=—  (35) 

L 

As  this  ratio  (SF)  becomes  1  or  smaller,  the  potential  of  unsatisfactory  perform¬ 
ance  in  terms  of  sliding  stability  of  a  structure  becomes  more  imminent. 
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The  unsatisfactory  performance  in  sliding  stability  for  a  structure  occurs 
along  a  theoretical  failure  surface,  which  is  identified  only  for  analysis  purposes 
and  does  not  necessarily  reflect  a  catastrophic  failure  of  the  structure.  The  failure 
surface  may  be  irregular  depending  on  the  homogeneity  of  the  backfill  and  the 
foundation  materials.  Therefore,  the  failure  surface  can  be  a  combination  of  fail¬ 
ure  planes  and  curved  surfaces.  The  stability  analysis  can  accordingly  be  reduced 
to  static  evaluation  of  assumed  wedges  of  failed  materials  with  assumed  force 
directions  at  the  interaction  surfaces  between  the  edges.  Therefore,  a  wedge 
configuration  is  determined  that  meets  the  conditions  of  statics.  Once  the  wedge 
configuration  is  determined,  the  R  and  L  terms  in  Equation  35  are  determined 
using  statics.  Then,  the  safety  factor  is  evaluated. 

The  analysis  in  CSLIDE  is  based  on  a  set  of  assumptions  (Pace  and  Noddin 
1987)  that  includes  (a)  wedge  interfaces  as  vertical  planes,  (b)  a  kinematically 
possible  failure  surface  that  consists  of  linear  segments,  (c)  two-dimensional  static 
analysis  of  edges,  (d)  the  shear  resistance  that  is  determined  using  Mohr-Coulomb 
failure  criteria,  (e)  force  equilibrium  that  is  satisfied,  (Q  moment  equilibrium  not 
considered,  (g)  negligible  shear  forces  at  wedge  interfaces,  (h)  wedges  that  have 
the  same  safety  factor,  (i)  small  displacements  that  do  not  have  any  effect  on 
resulting  wedge  forces,  and  (j)  the  stmcture  that  produces  one  structure  wedge 
because  of  its  relatively  high  shear  strength.  Usually,  an  iterative  computational 
procedure  is  needed  to  determine  the  SF  that  produces  force  equilibrium  for  the 
edges.  The  iterative  procedures  are  described  by  Pace  and  Noddin  (1987). 


CSLIDE  Program 

The  sliding  stability  of  concrete  structures  (CSLIDE)  program  (Pace  and 
Noddin  1987)  is  based  on  the  theory  discussed  in  the  preceding  section.  The 
program  was  written  in  FORTRAN  and  was  revised  to  run  in  a  batch  mode  using 
an  input  file.  An  output  file  of  CSLIDE  includes  the  safety  factor  as  defined  in 
Equation  35.  The  program  performs  determioistic  evaluation  of  SF  based  on  its 
input. 

An  input  file  of  CSLIDE  consists  of  the  following  information:  (a)  structural 
edge;  (b)  soil  on  both  sides  of  a  structure;  (c)  soil  below  the  structure,  (d)  methods 
of  analysis  that  include  single  plane,  multi-plane,  and  variable  angles;  (e)  water 
information;  (f)  input  wedge  angles;  and  (g)  loads.  Therefore,  CSLIDE  can 
handle  a  variety  of  soil  geometries,  structural  geometries,  and  loading  conditions. 

CSLIDE  has  the  following  limitations:  (a)  the  number  of  wedges  must  equal 
the  number  of  soil  layers;  (b)  soil  profiles  cannot  intersect;  (c)  once  the  base  of  a 
wedge  is  determined  in  a  soil  layer,  the  remaining  soil  profile  beyond  the  base  of 
the  wedge  is  ignored;  and  (d)  failure  occurs  from  left  to  right. 
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4  Reliability  Assessment 
for  Concrete  Gravity 
Structures 


The  reliability  of  sliding  stability  for  concrete  retaining  walls  and  gravity 
structures  can  be  assessed  using  the  reliability  methods  described  in  Chapter  2. 

In  the  first  section  of  this  chapter,  a  rehability  assessment  model  based  on  the 
ASM  method  is  introduced.  Numerical  algorithms  for  computing  derivatives  and 
finding  roots  of  equations  are  also  described.  In  the  second  section,  a  model  for 
nstng  DR  MCS  and  IS  is  described.  In  the  last  section,  the  development  of  the 
reliability  assessment  software  (RCSLIDE)  is  also  described. 


Performance  Function 

The  SF  for  the  sliding  stability  of  concrete  retaining  walls  and  gravity  stmc- 
tures  can  be  obtained  by  using  CSLIDE  (Pace  and  Noddin  1987)  and  can  be 
expressed  as 

SF=-  (36) 

L 


where 

R  =  resistance  in  the  form  of  a  restoring  strength 

L  =  load  effect  in  the  form  of  the  sliding  force 

From  Equation  36,  the  resistance  R  =  (SF)L .  Therefore,  the  performance 
function  Z  of  Equation  1  is  given  by 

Z=R-L=(SF)L-L=L(SF -1)  (37) 

Dividing  both  side  by  Z-  produces  the  following  revised  performance  function  G 
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G=  —  =  SF-1 
L 


(38) 


Therefore,  by  utilizing  CSLIDE  to  obtain  SF,  Equation  38  becomes  a  nonclosed 
expression  of  the  performance  function  as  follows: 

G=SF{X^,X2,-,X„)-1  (39) 

where  Z,  =  a  basic  strength  random  variable  such  as  geometry,  material  property, 
and  loads  that  are  input  quantities  to  CSLIDE  as  described  in  Chapter  3. 


Advanced  Second  Moment  (ASM)  Method 

The  ASM  method  as  described  in  Chapter  2  requires  a  closed-form  expression 
of  a  performance  function  (Equation  1).  In  Equation  39,  SF  is  not  in  a  closed- 
form  expression.  The  objective  of  this  section  is  to  introduce  the  needed  changes 
to  the  ASM  method  in  order  to  handle  nonclosed  forms  of  the  performance 
function. 

The  first  change  to  the  method  is  needed  in  computing  the  directional  cosines 
as  given  in  Equations  9,  21,  and  22  for  each  random  variable.  These  directional 
cosines  require  the  computation  of  the  partial  derivatives  of  G  with  respect  to  the 
random  variables.  The  partial  derivatives  with  respect  to  Xi  can  be  determined  as 

dG  _  dG  dSF  _(^.dSF  dSF  dSF(Xi,X2,-X„) 

dXi  dSF  dXi  dXi  dXi  dXi  ^ 

The  partial  derivative  with  respect  to  each  Xi  in  Equation  40  caimot  be  determined 
analytically.  Therefore,  a  numerical  method  for  computing  these  partial  deriva¬ 
tives  is  needed  and  is  described  later  in  this  chapter  as  Algorithm  4-1.  The  direc¬ 
tional  cosine  for  each  random  variable  can  now  be  determined  as  given  by 
Equations  9,  21,  and  22. 

The  second  change  to  the  method  is  needed  in  solving  for  the  reliability  index 
P  by  substituting  Equations  10,  23,  and  24  into  Equation  39.  This  solution 
requires  finding  the  root  for  the  following  expression: 


G  =  5F(Zi,Z2,-,X„)-1  =  G(P)  =  0  (41) 

where  G(P)  =  G  as  a  function  of  p.  Therefore,  a  numerical  method  for  finding  the 
roots  of  nonlinear  equations  is  needed  and  is  described  later  in  this  chapter  as 
Algorithms  4-2  and  4-3.  After  p  is  found,  the  design  point  is  updated  according  to 
Equations  10,  23,  and  24,  and  this  completes  Steps  5  and  6  in  Algorithms  2-1  and 
2-2. 
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Numerical  differentiation 


By  definition,  if  a  function  f(x)  is  continuous  in  the  domain  of  x,  then  the 
derivative  of  f(x)  can  be  expressed  as  follows: 

■i/W  (4; 

dx  Ax 

The  accuracy  of  computing  the  derivative  in  Equation  42  can  be  improved  as 
follows: 


df(x)  _  f(x  +  ^x)-  f(x-Ax) 
dx  2Ax 


(43) 


Therefore,  the  partial  derivative  for  each  Xi  in  Equation  40  can  be  expressed 
as 


dSF 

dXi 


lim 


SF(Xi,X2 


AX->-0 


(44) 

X,-  +AX,-,X„)-SF(Xi,X2,-,Xi  -A3r,-,X„) 

2AXi 


Because  this  numerical  approach  can  be  applied  to  both  closed  and  nonclosed 
expressions,  it  meets  the  purposes  of  computer  usage  for  reliability  assessment. 

Algorithm  4-1.  Based  on  Equation  9,  the  numerical  algorithm  for  computing 
the  partial  derivative  for  a  random  variable  Xi  in  SF  can  be  summarized  as 
follows: 


a.  Step  1:  Compute  X,-  +  AZ,-  ,  in  which  AZ,-  is  a  specified  very  small 
quantity  (or  increment). 

b.  Step!'.  Compute  =SF(Zj,Z2,"-,Z,- +AZ, •,•••, Z„)  using 

CSLIDE. 

c.  Step  3:  Compute  X^  -  AX^ ,  in  which  AX^  is  a  specified  very  small 
quantity  (or  increment). 

d.  Step4:  Compute  SF2  =5F(Zi,Z2,-",Z,- -AZ, •,•••, Z„)  using 
CSLIDE. 

e.  Step  5:  Estimate  dSF  /  dX^  «  (SF^  -SF2)/  2AZ,- . 

f.  Step  6:  Repeat  Steps  1  to  5  for  eachZ; . 


End  of  Algorithm  4-1. 
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Numerical  solution  for  roots  of  equations 

Newton’s  method  for  root  finding  is  commonly  used  for  this  purpose,  but  this 
method  does  not  always  converge  to  the  tme  values  for  some  cases  (Ayyub  and 
McCuen  1996).  Besides,  Newton’s  method  requires  the  evaluation  of  derivative. 
The  bisection  method  always  converges,  but  the  speed  of  finding  roots  within  an 
acceptable  tolerance  for  a  nonlinear  equation  can  be  slow.  Therefore,  a  root¬ 
finding  method  called  linear  interpolation  method  or  regula  falsi  method  (Gerald 
and  Wheatley  1984)  is  adopted  herein.  This  method  always  converges,  and  its 
convergence  rate  is  somewhat  slower  than  Newton’s  method  and  faster  than  the 
bisection  method. 

Assuming  that  a  function  f(x)  is  continuous  in  [a,  b],  and  /(«)  and  f(b) 
are  of  opposite  signs,  the  algorithm  for  the  regula  falsi  method  can  be  stated  as 
follows  (Gerald  and  Wheatley  1984): 

Algorithm  4-2. 

Do  While  (  I  b-a  |  a  tolerance_l  or  |  f  ( c )  |  a  tolerance_2 ) 
c=b-f (b)*(b-a)/(f (b)-f (a) ) 

If  f(c)*f(a)<0  Then 
b=c 

Else 

a=c 
End  If 
End  Do 

End  of  Algorithm  4-2. 

In  cases  where  f(x)  has  a  significant  curvature  between  a  and  b,  the  conver¬ 
gence  speed  becomes  slow  (Gerald  and  Wheatley  1984).  To  improve  convergence 
speed  in  this  case,  a  modified  linear  interpolation  method  can  be  used.  According 
to  the  modified  method,  the  value  of  f{x)  at  the  unchanged  end  position  is 

replaced  with  f(x)l2  as  described  in  the  following  revised  algorithm  (Gerald  and 
Wheatley  1984): 

Algorithm  4-3. 

Fa=f(a);  Fb=f(b);  FF=f(a) 

Do  While  (  I  b-a  |  a  tolerance_l  or  |  Fc  |  a  tolerance_2 ) 
c=b-Fb* ( b-a ) / ( Fb-Fa ) 

Fc=f (c) 

If  Fc*Fa  <  0  Then 
b=c 
Fb=Fc 

If  Fc*FF>0  Then  Fa=Fa/2 

Else 

a=c 

Fa=Fc 

If  Fc*FF  >  0  Then  Fb=Fb/2 

End  If 
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FF=Fc 
End  Do 

End  of  Algorithm  4-3. 

The  implementation  of  computer  programs  in  ASM  for  finding  p  was  based 
on  the  modified  linear  interpolation  method.  After  p  is  found,  the  design  point  is 
updated  according  to  Equations  10,  23,  and  24,  and  this  completes  Steps  5  and  6 
in  Algorithms  2-1  and  2-2. 


Flowchart 

The  procedure  for  reliability  assessment  of  concrete  retaining  walls  and 
gravity  structures  using  the  ASM  method  according  to  Algorithms  2-1,  2-2,  4-1, 
and  4-3  is  shown  in  Figure  1  in  the  form  of  a  flowchart.  The  computer  program 
developed  based  on  this  flowchart  results  in  the  reliability  index  P,  unsatisfactory 
performance  probability  ,  directional  cosine  a^- ,  partial  safety  factors  for  the 
input  random  variables,  and  the  design  point. 


Monte  Carlo  Simulation  (MCS)  Methods 

The  performance  function  in  Equation  40  can  also  be  used  in  MCS  methods 
for  reliability  assessment  of  concrete  retaining  walls  and  gravity  structures.  For 
the  purpose  of  IS,  Equation  40  does  not  indicate  which  random  variables  are  on 
the  resistance  side  (R)  and  which  random  variables  are  on  the  load-effect  side  (L). 
Therefore,  a  sensitivity  test  for  each  random  variable  is  needed  to  allocate  each 
random  variable  to  either  a  resistance  or  loading  side.  Some  random  variables 
might  not  belong  to  the  L  side  nor  to  the  R  side;  therefore,  they  can  be  assigned  to 
a  neutral  side  for  use  in  the  importance  sampling  procedure. 

Algorithm  4-4.  The  sensitivity  test  can  be  performed  using  the  following 
proposed  algorithm: 

a.  Step  1:  Compute  SF  =SF(Xi,X2,---,Xi,---,X„)  using  CSLIDE. 

b.  Step  2:  Compute  X^  +  AX^ ,  in  which  AXi  is  a  specified  very  small 
quantity  (or  increment). 

c.  Step 3:  Compute  SF^  =SF(X^,X2,-">Xi  +  AXi,---,Xj^)  using 
CSLIDE. 

d.  Step  4:  DIFFi  =  SFi  -  SF. 

e.  Step  5:  Compute  Z;  -  AZ,-,  in  which  AZ,-  is  a  specified  very  small 
quantity  (or  increment). 
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Figure  1 .  Advanced  second  moment  method  flowchart 
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f.  Step6\  Compute  SF2  =SF{X-^,X2,---,Xi- using 
CSLIDE. 

g.  Step  7:  DIFF2  =  SF2  -  SF 

h.  Step  8:  JfDIFFi  =  DIFF2,  ThenZ;  is  on  neutral  side,  Else  If  DZFFi  < 
DIFF2  Then  Z  is  on  loading  side,  Else  Z  is  on  resistant  side. 

i.  Step  9:  Repeat  Steps  2  to  8  for  each  Z  • 

End  of  Algorithm  4-4. 

After  the  sensitivity  tests  are  completed  for  all  the  random  variables,  the 
shifting  direction  of  the  mean  value  of  each  random  variable  can  be  determined  for 
the  importance  sampling  method.  A  random  number  generator  (Press  et  al.  1992) 
is  then  used  to  produce  random  numbers  as  CDF  values  for  the  corresponding 
basic  random  variables.  The  inverse  transformation  method  (Law  and  Kelton 
1991)  according  to  the  shifted  mean  values  and  unchanged  standard  deviations  is 
used  to  obtain  the  corresponding  generated  values  for  the  basic  variables. 

Estimates  of  the  probability  of  unsatisfactory  performance  and  its  coeffi¬ 
cient  of  variation  can  be  computed  using  Equations  28  and  29,  respectively. 

The  IS  method  requires  that  the  distributions  of  some  of  the  basic  random 
variables  be  shifted  according  to  Algorithm  4-4  in  order  to  increase  the  chance  of 
obtaining  unsatisfactory  performance  in  the  simulation  runs.  As  a  measure  of  the 
significance  of  shifting  these  variables  on  the  performance  of  the  IS  method,  a 
shifted  safety  factor  SFs  based  on  a  CSLIDE  run  using  the  shifted  mean  values  of 
the  basic  random  variables  can  be  utilized.  The  random  variables  can  be  shifted 
using  different  levels  that  depend  on  their  estimated  effects  on  SF  based  on  the 
sensitivity  analysis  of  Algorithm  4-4.  The  random  variables  need  to  be  incremen¬ 
tally  shifted  until  some  target  level  for  the  safety  factor  SF,  is  achieved  that  can 
ensure  the  desired  IS  performance.  In  the  case  of  correlated  random  variables, 
these  variables  should  not  be  selected  for  defining  the  sampling  (or  importance) 
density  function  %.  The  target  safety  factor  should  be  smaller  than  the  first-order 
mean  safety  factor  SF^  based  on  the  mean  values  of  the  unshifted  random  vari¬ 
ables.  The  smallest  value  the  target  safety  factor  can  take  is  one,  which  corre¬ 
sponds  to  a  point  on  the  unsatisfactory  performance  surface.  A  user  of  this 
method  needs  to  specify  a  target  safety  factor  based  on  judgment.  The  larger  the 
target  SF  value,  the  smaller  the  efficiency  of  the  IS  method.  On  the  other  hand, 
the  smaller  the  target  SF  value,  the  larger  the  shift  of  the  probability  distributions 
of  the  shifted  basic  random  variables,  which  might  result  in  violating  some  limits 
or  constraints  that  are  internal  to  the  CSLIDE  program.  A  large  shift  in  some 
random  variables  can  result  in  warning  messages  by  CSLIDE  due  to  violating 
some  internal  limits  and  constraints  within  CSLIDE;  e.g.,  water  levels  on  the  two 
sides  of  the  structure  are  below  the  stmcture  and  are  not  at  the  same  elevation. 
Therefore,  the  safety  factor  target  should  be  selected  carefully;  the  larger  the  target 
safety  factor,  the  better  IS  performance  with  a  penalty  of  an  increased  number  of 
simulation  cycles  needed.  Using  the  design  point  as  a  shifting  mean  for  the  target 
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safety  factor  is  not  necessarily  a  good  choice  for  IS  with  a  nonclosed  performance 
function. 

The  procedure  for  reliability  assessment  of  concrete  retaining  walls  and  grav¬ 
ity  structures  using  importance  sampling  and  direct  MCS  is  shown  in  Figure  2  in 
the  form  of  a  flowchart.  The  computer  program  developed  based  on  this  flowchart 
results  in  the  estimated  unsatisfactory  performance  probability. 


Development  of  RCSLIDE  Software 

The  software  RCSLIDE  (reliability  and  stability  assessment  of  concrete 
structures)  provides  stability  analysis  and  reliability  assessment  of  concrete 
retaining  walls  and  gravity  structures  with  a  user-friendly  interface.  RCSLIDE 
contains  the  following  three  main  analysis  functions:  (a)  stability  analysis  for 
concrete  retaining  walls  and  gravity  structures  using  CSLIDE,  (b)  reliability 
assessment  using  ASM  method,  and  (c)  reliability  assessment  using  MCS  with  IS 
and  DR.  In  addition  to  these  three  main  analysis  functions,  the  user  interface  of 
the  software  for  RCSLIDE  allows  the  preparation  of  input  data  for  running  these 
three  main  analysis  functions,  viewing  the  results  of  the  analysis,  and  plotting  the 
simulation  results  and  structural  geometry.  The  analysis  functions  were  all  written 
in  FORTRAN,  and  the  interface  portion  was  written  in  Microsoft  Visual  BASIC. 
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5  Case  Studies 


Two  case  studies  were  developed  in  this  study  for  noncoixelated  and  correlated 
random  variables.  The  primary  objective  of  these  case  studies  was  to  examine  the 
computational  accuracy  of  the  reliability  methods.  Probabilistic  characteristics  of 
random  variables  and/or  reliability  results  are  not  necessarily  representative  of  real 
gravity  structures. 


Noncorrelated  Random  Variables 

An  example  with  11  random  variables  is  described  in  this  section  to  illustrate 
reliability  assessment  of  the  sliding  stability  of  gravity  concrete  structures  using 
the  developed  software  RCSLIDE.  The  ASM  method,  MCS,  and  IS  are  used  in 
this  example. 

A  typical  gravity  concrete  structure  (Figure  3)  was  selected  from  the  example 
cases  in  the  user's  manual  for  CSLIDE  (Pace  and  Noddin  1987).  The  selected 
structure  involves  soil  and  water  pressure  on  both  its  sides  and  vertical  surcharges. 
Eleven  random  variables  were  then  identified  and  selected  from  the  input  to 
CSLIDE.  The  random  variables  were  chosen  from  structural,  sod,  and  loading 
properties.  These  variables  and  their  statistical  characteristics  were  selected  for 
the  purpose  of  demonstration,  and  they  do  not  necessarily  reflect  real  variables  or 
real  statistics.  The  11  random  variables  are  shown  in  Table  1.  Table  2  shows  the 
final  unsatisfactory  performance  probabdity  assessment  results  based  on  the  ASM 
method.  Figures  4-6  show  the  unsatisfactory  performance  probability  assessment 
results  based  on  simulation  methods  with  different  SF,  values  and  different  seed 
values.  The  results  for  both  ASM  and  MCS  are  summarized  in  Table  3.  The 
detailed  input  data  and  results  for  CSLIDE  and  RCSLIDE  are  provided  in 
Appendix  A. 

Table  3  and  Figures  4-6  show  that  the  results  using  ASM,  IS,  and  direct  MCS 
are  quite  in  agreement  except  for  the  results  based  on  importance  sampling  with  an 
SF,  of  1.0.  The  disagreement  in  this  case  can  be  attributed  to  the  large  shift  of 
random  variables  according  to  the  importance  function  that  resulted  in  an  SF,  of 
1.0.  The  large  shift  in  this  case  resulted  in  warning  messages  by  CSLIDE  due  to 
violating  some  internal  limits  and  constraints  within  CSLIDE;  e.g.,  water  levels  on 
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Figure  3.  Retaining  wall 


Table  1 
Characterisi 

tics  of  Input  Random  Varia 

bias 

X 

(1) 

Name 

(2) 

Meaning 

(3) 

Mean 

(4) 

COV(Stddv) 

(5) 

Type 

(6) 

1 

GAMC 

Equivalent  structure  unit  weight,  kef 

0.15 

0.15 

Normal 

2 

PHIL 

Angle  of  internal  friction,  deg  (left-side 
soil) 

28.00 

0.30 

Normal 

3 

GAML 

Layer  moist  unit  weight,  kef  (left-side 
soil) 

0.12 

0.10 

Normal 

4 

PHIR 

Angle  of  internal  friction,  deg  (right- 
side) 

28.00 

0.30 

Normal 

5 

GAMR 

Layer  moist  unit  weight,  kef  (right-side 
soil) 

0.12 

0.10 

Normal 

6 

PHIC 

Angle  of  internal  friction,  deg 
(structural  wedge  soil) 

30.00 

0.30 

Normal 

7 

WLL 

Left-side  water  elevation,  ft 

5.00 

(0.10) 

Normal 

8 

WLR 

Right-side  water  elevation,  ft 

1.50 

JOJO) _ 

Normal 

9 

PLO 

Point  load  magnitude,  kips 

2.00 

0.15 

Lognormal 

10 

SMAG 

Strip  load  magnitude,  kips/ft 

0.80 

0.15 

Lognormal 

11 

GRAM 

Ramp  load  magnitude,  kIps/ft 

0.55 

0.15 

Lognormal 

Note:  The  progra 
by  0.028.  Tocony 
convert  kips  to  kN 

m  RCSLIDE  uses  non-Sl  units  of  measurement.  To  convert  cu  ftto  cu  m,  multiply 
/ert  deg  to  rad,  multiply  by  0.0175.  To  convert  ft  to  m,  multiply  by  0.305.  To 
multipiv  bv  4.448.  To  convert  kips/ft  to  kN/m,  multiply  by  14.58. 

the  two  sides  of  the  structure  are  below  the  structure  and  are  not  at  the  same  eleva¬ 
tion.  Therefore,  the  safety-factor  target  should  be  selected  carefully;  the  larger  the 
target  safety  factor  the  better  IS  performance  with  a  penalty  of  an  increased  num¬ 
ber  of  needed  simulation  cycles. 
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Table  2 

Reliability  Results  Using  Advanced  Second  Moment  (ASM) 
Method 


X 

(1) 

Equivalent 
Normal  mean 
(2) 

Equivalent 
Normal  stddv 
(3) 

Directional 

Cosine 

(4) 

Design  Point 
(5) 

Partial  Safety 
Factor 
(6) 

1 

0.150000D+00 

0.225000D-01 

0.7200000-01 

0.1463060+00 

0.9753740+00 

2 

0.280000D4-02 

0,8400000+01 

0.6906080+00 

0.1477230+02 

0.5275840+00 

3 

0.120000D+00 

0.1 200000-01 

-0.7502140-01 

0.1220530+00 

0.1017110+01 

4 

0.280000D+02 

0.8400000+01 

0.3361170+00 

0.1220530+00 

0.7700760+00 

5 

0.120000D+00 

0.1200000-01 

0.5558700-01 

0.1184790+00 

0.9873250+00 

6 

0.300000D+02 

0.9000000+01 

0.6046320+00 

0.1759190+02 

0.5863960+00 

7 

0.500000D+01 

0.1000000+01 

-0.1461740-01 

0.5033330+01 

0.1006670+01 

8 

0.150000D+01 

0.1000000+01 

-0.3362160-01 

0.1576660+01 

0.1051110+01 

9 

0.197787D+01 

0.2950140+00 

0.0000000+00 

0.1977870+01 

0.9889360+00 

10 

0.789299D+00 

0.1100340+00 

0.1288840+00 

0.7569620+00 

0.9462020+00 

11 

0.542957D+00 

0.8599630-01 

-0.1124700+00 

0.5650110+00 

0.1027290+01 

Safety  factor  SF=  0,1 82545D+01 
Reliability  Index  p  =  0.228020D+01 

Unsatisfactory  performance  probability  Pu=  0.11 2980D-01 


Table  3 

Reliability  Results  Using  Advanced  Second  Moment  (ASM)  Method 
and  Monte  Carlo  Simulation  (MCS) 

Reliability  Method 

Unsatisfactory 

Performance 

Probability  Pu 

Reliability  Index  P 

ASM 

0.1129800-01 

0.2280200+01 

1  MCS  (50,000  cycles)  II 

IS01  -  SFt  =1 .0  (seed  =  -362646) 

0.9218470-03 

0.3114320+01 

IS02-  SFt  =1 .0  (seed  =  -41 751 44) 

0.2988460-02 

0.2749040+01 

IS01  -  SFt  =1 .5  (seed  =  -362646) 

0.1231900-01 

0.2247030+01 

IS02-  SFt=^.5  (seed  =  -4175144) 

0.1035430-01 

0.2313260+01 

OR01- SFf  =1.82545  (seed  =  -362646) 

0.0960000-01 

0.2341620+01 

OR02-  SF  =1.82545  (seed  =  -4175144) 

0.1000000-01 

0.2326350+01 

This  example  demonstrated  that  a  nonclosed  performance  function  can  be  used 
in  the  ASM  and  DR  and  IS  MCS  reliability  methods.  The  example  also  shows  the 
effect  of  selecting  a  seed  value  on  the  assessed  unsatisfactory  performance 
probability. 
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Figure  4.  IS  with  SFt  =1.0 


Correlated  Random  Variables 

This  case  study  documents  the  inclusion  of  correlated  random  variables  in 
RCSLIDE  and  their  effects  on  reliability  results.  This  second  example  uses  the 
same  retaining  wall  structure  as  shown  in  Figure  1;  however,  different  random 
variables  are  invoked  to  fully  show  the  effects  of  correlated  random  variables  in 
RCSLIDE. 

This  case  study  incorporates  negative  correlation  coefficients  between  the 
angle  of  internal  friction  (j)  and  cohesion  c  values  for  the  right,  left,  and  structural 
soil  layers.  The  correlation  coefficient  between  (j)  and  c  for  this  example  was 
determined  to  be  -0.5  and  the  distribution  of  ())  and  c  were  considered  normally 
distributed.  Note  that  RCSLIDE  requires  that  distributions  for  correlated  random 
variables  must  be  normally  distributed.  This  is  the  default  in  the  RCSLIDE 
program  and  the  user  is  alerted  to  this  requirement  prior  to  execution  of  the 
reliability  program.  The  characteristics  of  input  random  variables  are  shown  in 
Table  4.  The  input  files  for  both  the  CSLIDE  and  RCSLIDE  are  shown  in 
Appendix  B. 
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Figure  5.  IS  with  SFt  =  1 .50 


Results  of  the  ASM  method  for  this  case  study  with  and  without  correlation 
are  shown  in  Tables  5  and  6,  respectively.  The  results  show  that  the  use  of 
negative  correlation  for  soil  properties,  (j)  and  c,  resulted  in  a  slight  change  in  the 
reliability  index  for  the  gravity  structure  compared  with  the  reliability  index  for 
noncorrelated  soil  properties. 

The  results  from  the  direct  MCS  and  IS  compared  with  the  correlated  ASM 
are  shown  in  Table  7.  The  simulation  results  (direct  MCS  and  IS)  are  in  good 
agreement  with  the  results  of  the  ASM  method.  The  difference  between  simulation 
results  (direct  MCS  and  IS)  and  ASM  results  is  greater  for  direct  MCS  than  for  IS. 
This  can  be  directly  attributed  to  the  relatively  small  sample  size  (i.e.,  simulation 
cycles)  chosen  for  this  example.  The  complete  output  from  ASM  (correlated  and 
noncorrelated),  direct  MCS,  and  IS  is  shown  in  Appendix  B. 
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Figure  6.  Direct  MCS  with  SFt  =  1 .82545 


Table  4 

Characteristics  of  Input  Random  Variables  With  Correlation 


Meaning 


Equivalent  structure  unit  weight,  kef 


Angle  of  internal  friction,  deg  (left-side  soil 


Cohesion,  ksf  (left-side 


Laver  moist  unit  weight,  1 ,000  cu  ft  (left-side  soil 


Angle  of  Internal  friction,  deg  (right-side 


Cohesion,  ksf  (right-side 


Laver  moist  unit  weight,  1 .000  cu  ft  (right-side  soil 


Angle  of  internal  friction,  deg  (structural  wedge  soil 


Cohesion,  ksf  (structural  wedge 


Left-side  water  elevation,  ft  _ 


11  I  WLR  I  Right-side  water  elevation,  ft  1.50  (0.10)  Normal 

Note:  The  program  RCSLIDE  uses  non-SI  units  of  measurement.  To  convert  cu  ft  to  cu  m,  multiply  by  0.028.  To  convert  deg  to  rad, 
multiply  by  0.01 75.  To  convert  ft  to  m,  multiply  by  0.305.  To  convert  kips  to  kN,  multiply  by  4.448.  To  convert  kips/ft  to  kN/m,  multiply 
by  1 4.58.  -  - 
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Table  5 

Reliability  Results  Using  Advanced  Second  Moment  (ASM) 
Method  With  Correlation 


X 

(1) 

Equivalent 
Normal  mean 
(2) 

Equivalent 
Normal  stddv 
(3) 

Directional 

Cosine 

(4) 

Design  Point 
(5) 

Partial  Safety 
Factor 
(6) 

1 

0.150000D+00 

0.2250000-01 

0.7326000-01 

0.1454380+00 

0.9695840+00 

2 

0.280000D+02 

0.8400000+01 

0.2981470+00 

0.1505480+02 

0.5376730+00 

3 

0.500000D-01 

0.1000000-01 

0.470781  D+OO 

0.5715860-01 

0.1143170+01 

4 

0.120000D+00 

0.1200000-01 

-0.2028870+00 

0.1267390+00 

0.1056160+01 

5 

0,2800000+02 

0.8400000+01 

0.1676320+00 

0.2020510+02 

0.7216120+00 

6 

0.500000D-01 

0.1000000-01 

0.2903460+00 

0.5463980-01 

0.1092800+01 

7 

0.1 200000+00 

0.1200000-01 

0.7686910-01 

0.1174470+00 

0.9787240+00 

8 

0.3000000+02 

0.9000000+01 

0.3620250-00 

0.1196340+02 

0.3987810+00 

9 

0.5000000-01 

0.5000000-03 

0.6270460+00 

0.5050100-01 

0.1010020+01 

10 

0.5000000+01 

0.1000000+01 

-0.1960540-01 

0.5054260+01 

0.1010850+01 

11 

0.1000000+01 

0.1000000+01 

-0,0000000+00 

0.1000000+01 

0.1000000+01 

Reliability  index  P  =  0.276785D+01 

Unsatisfactory  performance  probability  Pu  =  0.2821 36D-02 


Table  6 

Reliability  Results  Using  Advanced  Second  Moment  (ASM) 
Method  Without  Correlation 


X 

(1) 

Equivalent 
Normal  mean 
(2) 

Equivalent 
Normal  stddv 
(3) 

Directional 

Cosine 

(4) 

Design  Point 
(5) 

Partial  Safety 
Factor 
(6) 

1 

0.1500000+00 

0.2250000-01 

0.5010940-01 

0.1469780+00 

0.9798560+00 

2 

0.2800000+02 

0.8400000+01 

0.6444040+00 

0.1349310+02 

0.4818950+00 

3 

0.5000000-01 

0.1000000-01 

0.3645870-01 

0.4902290-01 

0.9804580+00 

4 

0.1200000+00 

0.1200000-01 

-0.1674580+00 

0.1253850+00 

0.1044880+01 

5 

0.2800000+02 

0.8400000+01 

0.2834560+00 

0.2161880+02 

0.7721000+00 

6 

0.5000000-01 

0.1000000-01 

0.2793440-01 

0.4925140-01 

0.9850270+00 

7 

0.1200000+00 

0.1200000-01 

0.9240950-01 

0.1170280+00 

0.9752340+00 

8 

0.3000000+02 

0.9000000+01 

0.6744100+00 

0.1373310+02 

0.4577700+00 

9 

0.5000000-01 

0.5000000-03 

0.1363660-02 

0.4999820-01 

0.9999630+00 

10 

0.5000000+01 

0.1000000+01 

-0.2884841-01 

0.5077410+01 

0.1015480+01 

11 

0.1000000+01 

0.1000000+01 

-0.8676510-01 

0.1232530+01 

0.1232530+01 

Reliability  index  p  =  0.268002D+01 

Unsatisfactory  performance  probability  Pu  =  0.368090D-02 
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Table  7 

Reliability  Results  Using  Advanced  Second  Moment  (ASM)  Method 
and  Monte  Carlo  Simulation  (MCS)  With  Correlation 


Reiiabilitv  Method 

Unsatisfactory 

Perlormance 

Probability  Pu 

■■■ 

ASM 

0.2821 36D-02 

0.276785D+01  I 

MCS  fS.OOO  cycles)  _ _ _ _ 

IS-  SF,  =1.8  (seed  =  -9527654) 

0.28631 8D-02 

0.276305E+01 

DR- SF,  =1.931 98  (seed  =  -1142561) 

0.500000D-02 

0.257583E+01 
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6  Conclusions 


A  methodology  based  on  the  safety  factor  of  CSLIDE  was  developed  to  assess 
the  structural  reliability  for  concrete  retaining  walls  and  gravity  structures.  This 
study  demonstrates  that  the  structural  reliability  of  concrete  retaining  walls  and 
gravity  structures  can  be  assessed  with  a  nonclosed  performance  function  using 
the  advanced  second  moment  (ASM)  method,  direct  (DR)  Monte  Carlo  simulation 
(MCS),  and  importance  sampling  (IS).  Correlated  and  noncorrelated  random 
variables  were  considered  in  this  study. 

The  importance  sampling  simulation  requires  a  user-defined  target  safety 
factor  that  is  used  to  shift  the  mean  values  of  the  basic  random  variables.  In  this 
case,  the  design  point  is  not  necessarily  a  good  choice  for  shifting  the  mean  values, 
nor  is  a  target  safety  factor  of  one.  The  safety  factor  target  should  be  selected 
carefully;  the  larger  the  target  safety  factor,  the  better  IS  performance  with  a 
penalty  of  an  increased  number  of  simulation  cycles  needed.  Importance  sampling 
in  reliabihty  assessment  needs  to  be  carefully  employed  and  compared  with  other 
methods.  On  the  contrary,  direct  MCS  is  quite  good,  but  with  the  disadvantage  of 
possibly  a  large  number  of  simulation  cycles  in  order  to  achieve  an  accurate  esti¬ 
mated  unsatisfactory  performance  probability  for  structural  stability. 

Case  studies  were  used  to  demonstrate  utilizing  a  nonclosed  performance 
function  in  conjunction  with  the  ASM  and  DR  and  IS  MCS  reliability  methods. 
The  example  also  shows  the  effect  of  selecting  a  seed  value  and  a  target  safety 
factor  on  the  assessed  unsatisfactory  performance  probability.  In  addition,  the 
effect  of  correlation  of  random  variables  was  examined. 
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7  Recommended  Future 
Work 


In  this  study,  a  methodology  and  software  for  the  stability  and  reliability 
assessment  of  gravity  concrete  stmctures  were  developed.  The  software  devel¬ 
opment  procedure  described  in  this  study  is  significant  because  it  estabhshes  a 
prototype  reliability  software  that  is  modular  and  based  on  an  existing  CASE 
program.  Other  CASE  programs  can  be  modified  and  utilized  for  reliability 
purposes  in  a  similar  fashion.  The  following  recommendations  are  provided: 

a.  Case  studies  and  software  testing.  Case  studies  using  RCSLIDE  need  to 
be  developed.  Software  testing  can  also  be  performed. 

b.  Variable  arrays  and  their  values.  In  the  process  of  reliability  assessment, 
the  input  random  variables  used  by  CSLIDE  and  the  safety  factor  obtained 
from  CSLIDE  were  not  directly  passed  to  the  performance  function  in 
reliabihty  modules,  ASM  and  MCS.  Therefore,  speed  and  the  accuracy  in 
numerical  computations  were  somewhat  compromised.  Increased  speed 
and  accuracy  in  numerical  computations  without  the  loss  of  flexibility  for 
each  module  is  needed. 

c.  System  integration  and  help  menu.  The  CSLIDE  program  needs  to  be 
integrated  with  the  reliability  programs  resulting  in  one  integrated 
FORTRAN  program.  The  result  will  enhance  speed  and  performance. 

The  integrated  program  will  still  be  executed  from  the  Microsoft  Visual 
BASIC  shell. 

d.  Help.  To  reach  a  truly  user-friendly  interface,  a  hypertext-format-on-line 
help  manual  needs  to  be  developed  instead  of  the  current  text  format  help. 
The  help  capabilities  can  be  improved  to  allow  for  word  searching.  Also, 
the  help  file  can  be  converted  to  a  semi  or  full  hypertext  format  according 
to  Microsoft  Visual  BASIC  software. 

e.  Converters  or  ftlters  and  graphics.  Some  existing  input  data  files  for 
CSLIDE  use  free  format  and  cannot  be  used  directly  in  the  reliability 
assessment  computer  software  RCSLIDE,  which  uses  formatted  input  data 
as  shown  in  Appendix  A.  A  filter  program  for  converting  the  free  format 
data  to  formatted  data  is  needed  to  develop  and  facilitate  the  usage  of 
RCSLIDE.  A  filter  or  converter  can  be  developed  to  read  existing 
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CSLIDE  input  files  and  save  them  in  the  RCSLIDE  format.  The  filter  will 
read  an  existing  input  file  using  an  alphanumeric  format,  and  then  will 
strip  it  to  its  individual  components,  which  will  be  placed  at  the  correct 
locations  in  fields  of  RCSLIDE.  The  graphical  capabilities  of  RCSLIDE 
can  be  improved  to  include  plotting  results  from  its  output  file.  For 
example,  pressure  distributions  and  forces  can  be  added. 

f.  Soil  properties.  Modeling  correlated  soil  properties  should  be  improved 
by  developing  and  implementing  cases  involving  lognormal  distributions, 
with  and  without  normal  distributions,  for  soil  strength  parameters  such  as 
angle  of  internal  friction  and  cohesion  since  analytical  solution  methods 
exist  in  this  case.  Also,  the  nonlinear  relationships  for  soil  strength  prop¬ 
erties  should  be  developed  and  implemented  in  CSLIDE  instead  of  using 
soil  strength  parameters.  Then  RCSLIDE  can  be  revised  accordingly. 
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Appendix  A 

Input  and  Output  Data  for  Case 
Study  1 


Input  Data  for  CSLIDE 

START 

001  TITL  Title  Problem-l 
002  TITL 
003  TITL 
004  TITL 

005  STRU  8  1.500000E-01  O.OOOOOOE+00  l.OOOOOOE+00 

006  0.00  0.00 

007  0-00  2.00 

008  6.00  2-00 

009  6.00  14.00 

010  8-00  14.00 

011  8.00  2.00 

012  12-00  2.00 

013  12.00  0.00 

014  SOLT  1  1  2.800000E+01  O.OOOOOOE+00  1.200000E-01  1.400000E+01 

O.OOOOOOE+00 

015  -500.00  14.00 

016  SORT  1  1  2.800000E+01  O.OOOOOOE+00  1.200000E-01  4.000000E+00 

O.OOOOOOE+00 

017  500.00  4.00 

018  SOST  3-OOOOOOE+Ol  O.OOOOOOE+00  O.OOOOOOE+00 

019  METH  1 

020  WATR  5.000000E+00  1.500000E+00  6.250000E-02  -1 

021  FACT  0-5000  1.5000  1.0000 

022  VPLO  22.0000  2.000000E+00 

023  VSLO  10.0000  8.0000  8.000000E-01 

024  VRLO  6.0000  9.0000  5.500000E-01 

025  END 


Input  Data  for  Reliability  Methods 

TITL  Title  ASM-lbOl 

STRUOOOl  1.500000E-01  1.500000E-01  C  NOR 

SOLT0103  2.800000E+01  3.000000E-01  C  NOR 

SOLT0105  1.200000E-01  l.OOOOOOE-01  C  NOR 

SORT0103  2.800000E+01  3.000000E-01  C  NOR 

SORT0105  1.200000E-01  l.OOOOOOE-01  C  NOR 

SOST0003  3.000000E+01  3.000000E-01  C  NOR 

WATROOOl  5.000000E+00  l.OOOOOOE+00  S  NOR 

WATR0002  1.500000E+00  l.OOOOOOE+00  S  NOR 

VPL00105  2.000000E+00  1.500000E-01  C  LOG 


Appendix  A  Case  Study  1 


A1 


VSL00105  8.000000E-01  1.500000E-01  C  LOG 
VRLO0105  5.500000E~01  1.500000E-01  C  LOG 
END 


A2 


Output  of  Advanced  Second  Moment  (ASM) 
Method 


Title  =  Title  ASM-lbOl 

Total  number  of  random  variables  ~  11 

*****  The  input  data  ***** 


-x=; 

rv- 

— name- 

_ _ 

- mean - 

- Sigma - 

- cov - < 

code 

Type 

X( 

1) 

STRU 

0 

1 

.150000D+00 

.225000D-01 

.150000D+00 

C 

NOR 

X( 

2) 

SOLT 

1 

3 

.280000D+02 

.840000D+01 

.300000D+00 

C 

NOR 

X( 

3) 

SOLT 

1 

5 

.120000D+00 

.120000D-01 

.lOOOOOD+00 

C 

NOR 

X( 

4) 

SORT 

1 

3 

.280000D+02 

.840000D+01 

.300000D+00 

C 

NOR 

X( 

5) 

SORT 

1 

5 

.120000D+00 

-120000D-01 

.lOOOOOD+00 

C 

NOR 

X( 

6) 

SOST 

0 

3 

.300000D+02 

.900000D+01 

.300000D+00 

C 

NOR 

X 

7) 

WATR 

0 

1 

.500000D+01 

.lOOOOOD+01 

.200000D+00 

S 

NOR 

X 

8) 

WATR 

0 

2 

.150000D+01 

.lOOOOOD+01 

.666667D+00 

S 

NOR 

X( 

9) 

VPLO 

1 

5 

.200000D+01 

.300000D+00 

.150000D+00 

C 

LOG 

X( 

10) 

VSLO 

1 

5 

.800000D+00 

.120000D+00 

.150000D+00 

C 

LOG 

X( 

11) 

VRLO 

1 

5 

.550000D+00 

.825000D-01 

.150000D+00 

C 

LOG 

The 

tolerance 

of 

reliability 

index  (beta) 

=  .lOOOOOD 

-01 

The 

ratio  of  dx/x 

for  random 

variables 

=  .lOOOOOD 

-02 

The 

max.  iterations  for  finding  beta 

=  10 

The 

increment 

for 

finding  beta 

=  .2500 

The 

1= 

absolute  2= 

relative  tolerance  for  beta  =  2 

The 

limit  state  Z 

=R-L=L(SF-1) 

or  z==SF-l 

*****  The  result  ***** 


The  iteration  No.  =  1 


-i - meanEN - 

1  .150000D+00 

2  .280000D+02 

3  .120000D+00 

4  .280000D+02 

5  .120000D+00 

6  .300000D+02 

7  .500000D+01 

8  .150000D+01 

9  .197775D+01 

10  .791100D+00 

11  .543881D+00 
The  reliability 


- stdEN - 

.225000D-01 
.840000D+01 
.120000D-01 
.840000D+01 
.120000D-01 
.900000D+01 
.lOOOOOD+01 
.lOOOOOD+01 
.298333D+00 
.119333D+00 
.820415D-01 
index  ( beta ) 


— alpha 
.105081D+00 
.639212D+00 
•.119625D+00 
.329423D+00 
.671452D-01 
.611215D+00 
-.542979D-01 
•.896885D-01 
•OOOOOOD+OO 
.195616D+00 
•  .173740D+00 
.228933D+01 


- dp - 

.144587D+00 

.157077D+02 

.123286D+00 

.216651D+02 

.118155D+00 

.174066D+02 

.512431D+01 

.170533D+01 

.197775D+01 

.737659D+00 

.576513D+00 


-Partial. SF- 
.963915D+00 
.560990D+00 
.102739D+01 
.773753D+00 
.984628D+00 
.580219D+00 
.102486D+01 
.113688D+01 
.988875D+00 
.922074D+00 
.104821D+01 


The  iteration  No.  =  2 


-i  - 

— meanEN - 

- StdEN - 

1 

.150000D+00 

.225000D-01 

2 

.280000D+02 

.840000D+01 

3 

.120000D+00 

.120000D-01 

4 

.280000D+02 

.840000D+01 

5 

.120000D+00 

.120000D-01 

6 

.300000D+02 

.900000D+01 

7 

.500000D+01 

-lOOOOOD+01 

8 

.150000D+01 

.lOOOOOD+01 

9 

.197787D+01 

.295014D+00 

10 

.789299D+00 

.110034D+00 

11 

-542957D+00 

.859963D-01 

The 

reliability 

index  (beta) 

— alpha 
.720000D-01 
.690608D+00 
.750214D-01 
.336117D+00 
.555870D-01 
.604632D+00 
.146174D-01 
.336216D-01 
.OOOOOOD+OO 
.128884D+00 
.112470D+00 
.228020D+01 


- dp - 

.146306D+00 

.147723D+02 

.122053D+00 

.215621D+02 

.118479D+00 

.175919D+02 

.503333D+01 

.157666D+01 

.197787D+01 

.756962D+00 

.565011D+00 


-Partial. SF- 
.975374D+00 
.527584D+00 
.101711D+01 
.770076D+00 
.987325D+00 
.586396D+00 
.100667D+01 
.105111D+01 
.988936D+00 
.946202D+00 
.102729D+01 


The  unsatisfactory  performance  probability  =  .112980D-01 
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Output  of  Importance  Sampling  -  IS01  {SFt  =1 .0) 

Title  =  Title  ASM-lbOl 

Total  number  of  random  variables  =  11 

*****  The  input  data  ***** 


-x= 

=rv- 

— name- 

— 

- mean - 

- Sigma-^ — 

- cov - 

code 

Type 

X( 

1) 

STRU 

0 

1 

.150000D+00 

.2250000-01 

.1500000+00 

C 

NOR 

X( 

2) 

SOLT 

1 

3 

,2800000+02 

.8400000+01 

.3000000+00 

C 

NOR 

X( 

3) 

SOLT 

1 

5 

.1200000+00 

,1200000-01 

.1000000+00 

C 

NOR 

X( 

4) 

SORT 

1 

3 

.280000D+02 

.8400000+01 

,3000000+00 

C 

NOR 

X( 

5) 

SORT 

1 

5 

.1200000+00 

,1200000-01 

.1000000+00 

C 

NOR 

X( 

6) 

SOST 

0 

3 

.3000000+02 

.9000000+01 

,3000000+00 

C 

NOR 

X( 

7) 

WATR 

0 

1 

.5000000+01 

.1000000+01 

.2000000+00 

S 

NOR 

X( 

8) 

WATR 

0 

2 

,1500000+01 

.1000000+01 

,6666670+00 

S 

NOR 

X( 

9) 

VPLO 

1 

5 

.2000000+01 

.3000000+00 

.1500000+00 

C 

LOG 

X( 

10) 

VSIiO 

1 

5 

,8000000+00 

.1200000+00 

.1500000+00 

C 

LOG 

X( 

11) 

VRLO 

1 

5 

.5500000+00 

.8250000-01 

.1500000+00 

C 

LOG 

Tolerance  for  root  finding  in  beta  =  -lOOOOOD-Ol 
Tolerance  for  finding  dz/dx  of  general  rv  =  ,1000000-02 
Max.  iterations  for  finding  safety  index  =  10 


The  limit  state  Z=R-L=L(SF-1) 

*****  The  result  ***** 

The  simulation  cycle  = 

50000 

The  initial  given  seed  = 

0 

The  sensitivity  test  ratio 

=  .1000000+00 

The  target  safety  factor 

=  .1000000+00 

The  initial  output  cycle  increment  =  500 

«<  The  sensitivity  result  »> 


The  safety  factor  (SF)  =  .182545E+01 


-x=rv-  — Name- 

Sense. — 

— 

— M=Mean — 

SF1:M*  .90 

SF2:M*1.10 

— SFl-SF-  - 

— SF2-SF— 

X(  1)  STRU 
.90000+00 

0 

1 

.15000+00 

.17960+01 

.18540+01 

-.28960-01 

.28860-01 

X(  2)  SOLT 
.90000+00 

1 

3 

.28000+02 

.17390+01 

.19150+01 

-.86800-01 

.89260-01 

X{  3)  SOLT 
.11000+01 

1 

5 

.12000+00 

.18780+01 

.17790+01 

.52980-01 

-.46150-01 

X(  4)  SORT 
.90000+00 

1 

3 

.28000+02 

.17820+01 

.18720+01 

-.43900-01 

.46840-01 

X(  5)  SORT 
.90000+00 

1 

5 

.12000+00 

.17980+01 

.18530+01 

-.27580-01 

.27810-01 

X(  6)  SOST 
.90000+00 

0 

3 

.30000+02 

.17430+01 

.19120+01 

-.82060-01 

.86330-01 

X(  7)  WATR 
.11000+01 

0 

1 

.50000+01 

.18370+01 

.18140+01 

.11170-01 

-.11260-01 

X(  8)  WATR 
.11000+01 

0 

2 

.15000+01 

.18310+01 

.18200+01 

.55010-02 

-.56180-02 

X(  9)  VPLO 
.10000+01 

1 

5 

.20000+01 

.18250+01 

.18250+01 

.00000+00 

.00000+00 

X(  10)  VSLO 
.90000+00 

1 

5 

.80000+00 

.17720+01 

.18800+01 

-.53660-01 

.54580-01 

X(  11)  VRLO 
.11000+01 

1 

5 

.55000+00 

.18760+01 

.17790+01 

.50200-01 

-.46110-01 

«<  The  shifted  means  »> 


-x= 

=rv- 

— Name- 

— 

- M=Mean - 

Shifted. Mean 

- Ratio - 

X( 

1) 

STRU 

0 

1 

.1500000+00 

.1093500+00 

.7290000+00 

X( 

2) 

SOLT 

1 

3 

.2800000+02 

.2041200+02 

.7290000+00 

X( 

3) 

SOLT 

1 

5 

.1200000+00 

.1597200+00 

.1331000+01 

X( 

4) 

SORT 

1 

3 

.2800000+02 

.2041200+02 

.7290000+00 

X( 

5) 

SORT 

1 

5 

.1200000+00 

.8748000-01 

.7290000+00 
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X{  6)  SOST  0  3 

X(  7)  WATR  0  1 

X(  8)  WATR  0  2 

X(  9)  VPLO  1  5 

X(  10)  VSLO  1  5 

X(  11)  VRLO  1  5 


.300000D+02 

.500000D+01 

.150000D+01 

.200000D+01 

.800000D+00 

.550000D+00 


.218700D+02 

.665500D+01 

.199650D+01 

.200000D+01 

.583200D+00 

.732050D+00 


.729000D+00 

.133100D+01 

.133100D+01 

.lOOOOOD+01 

.729000D+00 

.133100D+01 


The  safety  factor  =  .930955D+00 

The  iteration  of  shifting  -  3 

The  initial  random  seed  =  -362646 


***  The  result  of  importance  sampling  simulation  *** 


No. of. cycle-  -ave.cum.Pu- 


500 

.100956E-03 

1000 

.196125E-03 

1500 

.456397E-03 

2000 

.396067E-03 

2500 

.322244E-03 

3000 

.274498E-03 

3500 

.292017E-03 

4000 

.289559E-03 

4500 

.266437E-03 

5000 

.241174E-03 

5500 

.242164E-03 

6000 

.295676E-03 

6500 

.309899E-03 

7000 

,289754E-03 

7500 

.279749E-03 

8000 

.275022E-03 

8500 

.259333E-03 

9000 

.245442E-03 

9500 

.234040E-03 

10000 

.223853E-03 

10500 

.214773E-03 

11000 

.207800E-03 

11500 

.293170E-03 

12000 

.296667E-03 

12500 

.288635E-03 

13000 

.350232E-03 

13500 

.339149E-03 

14000 

.334743E-03 

14500 

.342194E-03 

15000 

.331299E-03 

15500 

.321954E-03 

16000 

.315085E-03 

16500 

.310017E-03 

17000 

.311496E-03 

17500 

.303805E-03 

18000 

.295958E-03 

18500 

.295868E-03 

19000 

.288948E-03 

19500 

.281823E-03 

20000 

.280597E-03 

20500 

.275871E-03 

21000 

.269890E-03 

21500 

.290372E-03 

22000 

.147684E-02 

22500 

.144566E-02 

23000 

.141676E-02 

23500 

.139081E-02 

24000 

.136250E-02 

24500 

.133629E-02 

25000 

.130982E-02 

25500 

.129265E-02 

26000 

.127804E-02 

26500 

.125415E-02 

27000 

.123323E-02 

27500 

.121099E-02 

28000 

.118997E-02 

28500 

.123509E-02 
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29000  .121432E-02 
29500  .119612E-02 
30000  .117662E-02 
30500  .115807E~02 
31000  .114949E-02 
31500  .113138E-02 
32000  .111411E-02 
32500  .125021E-02 
33000  .123166E-02 
33500  .121409E-02 
34000  .119670E-02 
34500  .117969E-02 
35000  .116477E-02 
35500  .114886E-02 
36000  .113348E-02 
36500  .111863E-02 
37000  .110425E-02 
37500  .108955E-02 
38000  .110085E-02 
38500  .108839E-02 
39000  .107457E-02 
39500  .106399E-02 
40000  .105695E-02 
40500  .104511E-02 
41000  .104004E-02 
41500  .105100E-02 
42000  .103884E-02 
42500  .102745E-02 
43000  .104360E-02 
43500  .103304E-02 
44000  .102611E-02 
44500  .101544E-02 
45000  .100652E-02 
45500  .995545E-03 
46000  .990373E-03 
46500  .980915E-03 
47000  .971398E-03 
47500  .961226E-03 
48000  .952945E-03 
48500  .943868E-03 
49000  .934695E-03 
49500  .925331E-03 
50000  .921847E-03 


The  importance  failure  count  =  28770 
The  safety  factor  >  10  count  =  277 
The  reliability  index  =  .3114320+01 
The  binomial  statitics 

The  ave  of  Pu  (AvePu)  =  .9218470-03 
The  var  of  avePu  (varAvePu)  =  .1841990-07 
The  std  of  avePu  (stdAvePu)  =  .1357200-03 
The  cov  of  avePu  (covAvePu)  =  .1472260+00 


The  simulation  and  sample  mean  statistics 
The  ave  of  Pu  (avePu)  =  .9218470-03 

The  var  of  Pu  (varPu)b  =  .1442350-01 

The  var  of  Pu  (varPu)  =  .1442380-01 

The  std  of  Pu  (stdPu)  =  .1200990+00 

The  cov  of  Pu  (covPu)  =  .1302810+03 

The  var  of  avePu  (var AvePu )b  =  .2884700-06 

The  var  of  avePu  (var AvePu)  =  .2884760-06 

The  std  of  avePu  (stdAvePu)  =  .5370990-03 

The  cov  of  avePu  (covAvePu)  =  .5826340+00 

b  =  biased,  for  variance  computation  using  the  sample  size  N. 
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Output  of  Importance  Sampling  -  IS02  (SFf  =1.0) 

Title  =  Title  ASM-lbOl 

Total  number  of  random  variables  =  11 


*****  The  input  data  ***** 


X( 

X( 

X( 

X( 

X( 

X( 

X( 

X( 


:rv- 

— name- 

_ 

- mean - 

- Sigma - 

- cov - 

code 

Type 

1) 

STRU 

0 

1 

.1500000+00 

.2250000-01 

.1500000+00 

C 

NOR 

2) 

SOLT 

1 

3 

.2800000+02 

-8400000+01 

.3000000+00 

C 

NOR 

3) 

SOLT 

1 

5 

.1200000+00 

.1200000-01 

.lOOOOOO+OO 

C 

NOR 

4) 

SORT 

1 

3 

.2800000+02 

.8400000+01 

.300000O+00 

C 

NOR 

5) 

SORT 

1 

5 

.1200000+00 

.120000O-01 

.lOOOOOO+OO 

C 

NOR 

6) 

SOST 

0 

3 

.3000000+02 

.9000000+01 

.3000000+00 

C 

NOR 

7) 

WATR 

0 

1 

.5000000+01 

.1000000+01 

.2000000+00 

S 

NOR 

8) 

WATR 

0 

2 

.1500000+01 

.1000000+01 

.6666670+00 

S 

NOR 

9) 

VPLO 

1 

5 

.2000000+01 

.3000000+00 

.1500000+00 

C 

LOG 

10) 

VSLO 

1 

5 

.8000000+00 

.1200000+00 

.1500000+00 

C 

LOG 

11) 

VRLO 

1 

5 

.5500000+00 

.8250000-01 

.1500000+00 

C 

LOG 

The  limit  state  Z=R-“L=L(SF--1 ) 


*****  The  result  ***** 


The  simulation  cycle  =  50000 
The  initial  given  seed  =  0 
The  sensitivity  test  ratio  = 
The  target  safety  factor  = 
The  initial  output  cycle  increment  = 


.lOOOOOD+00 

.lOOOOOD+00 

500 


«<  The  sensitivity  result  »> 


The  safety  factor 

(SF)  -  . 

182545E+01 

-x=rv-  — Name- 

_ _ 

-M=Mean — 

SFlzM*  .90 

SF2;M*1.10 

— SFl-SF — 

— SF2-SF — 

Sense. — 

X(  1)  STRU 

0 

1 

.15000+00 

.17960+01 

.18540+01 

-.28960-01 

.28860-01 

.90000+00 

X(  2)  SOLT 

1 

3 

.28000+02 

.17390+01 

.19150+01 

-.86800-01 

.89260-01 

.90000+00 

X(  3)  SOLT 

1 

5 

.12000+00 

.18780+01 

.17790+01 

.52980-01 

-.46150-01 

.11000+01 

X(  4)  SORT 

1 

3 

.28000+02 

.17820+01 

.18720+01 

-.43900-01 

.46840-01 

.90000+00 

X(  5)  SORT 

1 

5 

.12000+00 

.17980+01 

.18530+01 

-.27580-01 

.27810-01 

.90000+00 

X(  6)  SOST 

0 

3 

.30000+02 

.17430+01 

.19120+01 

-.82060-01 

.86330-01 

.90000+00 

X(  7)  WATR 

0 

1 

.50000+01 

.18370+01 

.18140+01 

.11170-01 

-.11260-01 

.11000+01 

X(  8)  WATR 

0 

2 

.15000+01 

.18310+01 

.18200+01 

.55010-02 

-.56180-02 

.11000+01 

X(  9)  VPLO 

1 

5 

.20000+01 

.18250+01 

.18250+01 

.00000+00 

.00000+00 

.10000+01 

X{  10)  VSLO 

1 

5 

.80000+00 

.17720+01 

.18800+01 

-.53660-01 

.54580-01 

.90000+00 

X{  11)  VRLO 

1 

5 

.55000+00 

.18760+01 

.17790+01 

.50200-01 

-.46110-01 

.11000+01 


«<  The  shifted  means  »> 

-x=rV“  — Name — - M=Mean - 

X(  1)  STRU  0  1  .1500000+00 

X(  2)  SOLT  1  3  .280000O+02 

X(  3)  SOLT  1  5  .120000D+00 

X(  4)  SORT  1  3  .2800000+02 

X(  5)  SORT  1  5  .1200000+00 

X(  6)  SOST  0  3  .300000O+02 

X(  7)  WATR  0  1  .5000000+01 

X(  8)  WATR  0  2  .1500000+01 

X(  9)  VPLO  1  5  .2000000+01 


Shifted. Mean 
.1093500+00 
.2041200+02 
.1597200+00 
.2041200+02 
.8748000-01 
.2187000+02 
.6655000+01 
.1996500+01 
.2000000+01 


- Ratio - 

.7290000+00 

.7290000+00 

.1331000+01 

.7290000+00 

.729000D+00 

.7290000+00 

.1331000+01 

.1331000+01 

.1000000+01 
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X(  10)  VSLO  1  5  .800000D+00  .583200D+00  .729000D+00 
X(  11)  VRLO  1  5  .550000D+00  .732050D+00  .133100D+01 

The  safety  factor  =  .930955D+00 

The  iteration  of  shifting  =  3 

The  initial  random  seed  =  -4175144 

***  The  result  of  importance  sampling  simulation  *** 


No. of. cycle-  -ave.cum.Pu- 


500 

.119824E-03 

1000 

.629463E-04 

1500 

.261312E-03 

2000 

.197144E-03 

2500 

.168649E-03 

3000 

.173870E-03 

3500 

.157569E-03 

4000 

.236490E-03 

4500 

.227291E-03 

5000 

,209060E-03 

5500 

.224515E-03 

6000 

.206776E-03 

6500 

.227723E-03 

7000 

.259903E-03 

7500 

.301265E-03 

8000 

.285042E-03 

8500 

.269437E-03 

9000 

.387376E-03 

9500 

.376566E-03 

10000 

.385317E-03 

10500 

.367690E-03 

11000 

.357962E-03 

11500 

.346093E-03 

12000 

.334930E-03 

12500 

.325378E-03 

13000 

.322924E-03 

13500 

.313199E-03 

14000 

.302524E-03 

14500 

.292869E-03 

15000 

-283279E-03 

15500 

.274190E-03 

16000 

.268454E-03 

16500 

.620358E-03 

17000 

.615184E-03 

17500 

.605702E-03 

18000 

.592766E-03 

18500 

.607881E-03 

19000 

.639766E-03 

19500 

-563074E-02 

20000 

.549124E-02 

20500 

.536162E-02 

21000 

.523449E-02 

21500 

.511659E-02 

22000 

.500623E-02 

22500 

.496222E-02 

23000 

.485512E-02 

23500 

.475240E-02 

24000 

.465357E-02 

24500 

.456086E-02 

25000 

.448282E-02 

25500 

.439813E-02 

26000 

.431716E-02 

26500 

.423777E-02 

27000 

.416161E-02 

27500 

.408825E-02 

28000 

.401980E-02 

28500 

.394970E-02 

29000 

.388260E-02 

29500 

.381746E-02 

30000 

.375398E-02 

30500 

.369433E-02 
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31000  .437118E-02 
31500  .432964E-02 
32000  .426240E-02 
32500  .420035E-02 
33000  .414427E-02 
33500  .408270E-02 
34000  .403066E-02 
34500  .397249E-02 
35000  .392409E-02 
35500  .386926E-02 
36000  .381591E-02 
36500  .376435E-02 
37000  .371404E-02 
37500  .366497E-02 
38000  .362055E-02 
38500  .357532E-02 
39000  .372158E-02 
39500  .367979E-02 
40000  .363453E-02 
40500  .360308E-02 
41000  .356096E-02 
41500  .351836E-02 
42000  .347662E-02 
42500  .343614E-02 
43000  .339885E-02 
43500  .336030E-02 
44000  .332577E-02 
44500  .328947E-02 
45000  .330596E-02 
45500  .327190E-02 
46000  .323804E-02 
46500  .320370E-02 
47000  .316994E-02 
47500  .313851E^02 
48000  .310588E-02 
48500  .307394E-02 
49000  .304727E-02 
49500  .301671E-02 
50000  .298846E-02 


The 

importance  failure  count 

ss 

28881 

The 

safety  factor 

>  10  count 

= 

267 

The 

reliability  index 

- 

.274904D+01 

The 

binomial  statistics 

The 

ave  of  Pu 

(AvePu) 

-298846D-02 

The 

var  of  avePu 

(varAvePu) 

.595906D-07 

The 

std  of  avePu 

( stdAvePu ) 

= 

.244112D-03 

The 

cov  of  avePu 

(covAvePu) 

= 

.816848D-01 

The 

simulation  and  sample  mean 

statistics 

The 

ave 

of 

Pu 

( avePu ) 

= 

.298846D-02 

The 

var 

of 

Pu 

( varPu ) b 

.202745D+00 

The 

var 

of 

Pu 

(varPu) 

= 

.2027490+00 

The 

std 

of 

Pu 

( stdPu) 

= 

.4502760+00 

The 

cov 

of 

Pu 

( covPu ) 

= 

.1506720+03 

The 

var 

of 

avePu 

( varAvePu ) b 

= 

.4054900-05 

The 

var 

of 

avePu 

( varAvePu ) 

= 

.4054980-05 

The 

std 

of 

avePu 

( StdAvePu ) 

= 

.2013700-02 

The 

cov 

of 

avePu 

( covAvePu ) 

.6738240+00 

b  = 

biased, 

for  variance  computation  using 

the  sample  size  N. 
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Output  of  Importance  Sampling  -  IS01  {SFt  =1 .5) 

Title  =  Title  ASM-lbOl 

Total  number  of  random  variables  =  11 


*****  I 

-x=rv- 

The  input 

— name — 

data  ***** 

- mean - - 

- Sigma - 

- cov - 

code 

Type 

X( 

1) 

STRU 

0 

1 

.1500000+00 

.2250000-01 

.1500000+00 

C 

NOR 

X( 

2) 

SOLT 

1 

3 

.2800000+02 

.8400000+01 

.3000000+00 

C 

NOR 

X( 

3) 

SOLT 

1 

5 

.1200000+00 

.1200000-01 

.1000000+00 

C 

NOR 

X( 

4) 

SORT 

1 

3 

.2800000+02 

.8400000+01 

.3000000+00 

C 

NOR 

X( 

5) 

SORT 

1 

5 

.1200000+00 

.1200000-01 

.1000000+00 

C 

NOR 

X( 

6) 

SOST 

0 

3 

.3000000+02 

.9000000+01 

.3000000+00 

C 

NOR 

X( 

7) 

WATR 

0 

1 

.5000000+01 

.1000000+01 

.2000000+00 

S 

NOR 

X( 

8) 

WATR 

0 

2 

.1500000+01 

.1000000+01 

.6666670+00 

S 

NOR 

X( 

9) 

VPLO 

1 

5 

.2000000+01 

.300000D+00 

.1500000+00 

C 

LOG 

X( 

10) 

VSLO 

1 

5 

.8000000+00 

.1200000+00 

.1500000+00 

C 

LOG 

X( 

11) 

VRLO 

1 

5 

.5500000+00 

.8250000-01 

.1500000+00 

C 

LOG 

The 

1  limit  state 

Z=R-L=L(SF-1) 

or  Z=SF-1 

The 

The 

The 

The 

The 


simulation  cycle  =  50000 

initial  given  seed  =  -362646 

sensitivity  test  ratio  = 

target  safety  factor  = 

initial  output  cycle  increment  = 


lOOOOOD+00 

150000D+01 

0 


«<  The  sensitivity  result  »> 


The  safety  factor  (SF)  =  .182545E+01 


-x=rv-  — Name- 
Sense  .  — 

— 

— M=Mean — 

SF1:M*  .90 

SF2:M*1.10 

— SFl-SF-  - 

— SF2-SF— 

X(  1)  STRU 
.90000+00 

0 

1 

.15000+00 

.17960+01 

.18540+01 

-.28960-01 

.28860-01 

X(  2)  SOLT 
.90000+00 

1 

3 

.28000+02 

.17390+01 

.19150+01 

-.86800-01 

.89260-01 

X(  3)  SOLT 
.11000+01 

1 

5 

.12000+00 

.18780+01 

.17790+01 

.52980-01 

-.46150-01 

X(  4)  SORT 
.90000+00 

1 

3 

.28000+02 

.17820+01 

.18720+01 

-.43900-01 

.46840-01 

X(  5)  SORT 
.90000+00 

1 

5 

.12000+00 

.17980+01 

.18530+01 

-.27580-01 

.27810-01 

X(  6)  SOST 
.90000+00 

0 

3 

.30000+02 

.17430+01 

.19120+01 

-.82060-01 

.86330-01 

X(  7)  WATR 
.11000+01 

0 

1 

.50000+01 

.18370+01 

.18140+01 

.11170-01 

-.11260-01 

X(  8)  WATR 
.11000+01 

0 

2 

.15000+01 

.18310+01 

.18200+01 

.55010-02 

-.56180-02 

X(  9)  VPLO 
.10000+01 

1 

5 

.20000+01 

.18250+01 

.18250+01 

.00000+00 

.00000+00 

X(  10)  VSLO 
.90000+00 

1 

5 

.80000+00 

.17720+01 

.18800+01 

-.53660-01 

.54580-01 

X(  11)  VRLO 

1 

5 

.55000+00 

.18760+01 

.17790+01 

.50200-01 

-.46110-01 

.llOOD+01 


«<  The  shifted  means  »> 


-x= 

=rv- 

— Name- 

— 

- M=Mean - 

Shifted. Mean 

- Ratio - 

X( 

1) 

STRU 

0 

1 

.1500000+00 

.1350000+00 

.9000000+00 

X( 

2) 

SOLT 

1 

3 

.2800000+02 

.2520000+02 

.9000000+00 

X{ 

3) 

SOLT 

1 

5 

.1200000+00 

.1320000+00 

.1100000+01 

X( 

4) 

SORT 

1 

3 

.2800000+02 

.2520000+02 

.9000000+00 

X{ 

5) 

SORT 

1 

5 

.1200000+00 

.1080000+00 

.900000D+00 

X( 

6) 

SOST 

0 

3 

.3000000+02 

.2700000+02 

.9000000+00 

X( 

7) 

WATR 

0 

1 

.5000000+01 

.5500000+01 

.1100000+01 

X( 

8) 

WATR 

0 

2 

.1500000+01 

.1650000+01 

.1100000+01 

X( 

9) 

VPLO 

1 

5 

.2000000+01 

.2000000+01 

.1000000+01 

X( 

10) 

VSLO 

1 

5 

.8000000+00 

.7200000+00 

.900000D+00 

X( 

11) 

VRLO 

1 

5 

.5500000+00 

.6050000+00 

.1100000+01 

Appendix  A  Case  Study  1 


A9 


The  safety  factor  =  .143618D+01 

The  iteration  of  shifting  =  1 

****•>:  The  result  ***** 

The  revised  output  cycle  increment  =  500 

***  The  result  of  importance  sampling  simulation  *** 


No. of. cycle-  -ave.cum.Pu- 


500 

.949377E-02 

1000 

.823076E-02 

1500 

.127657E-01 

2000 

.105944E-01 

2500 

.963586E-02 

3000 

.977791E-02 

3500 

.106235E-01 

4000 

.103422E-01 

4500 

.977378E-02 

5000 

.979711E-02 

5500 

.966660E-02 

6000 

.113069E-01 

6500 

.124313E-01 

7000 

.118592E-01 

7500 

.115973E-01 

8000 

.121058E-01 

8500 

.118632E-01 

9000 

.120089E-01 

9500 

.116215E-01 

10000 

.120883E-01 

10500 

.120262E-01 

11000 

.118223E-01 

11500 

.119256E-01 

12000 

.121193E-01 

12500 

.118127E-01 

13000 

.116685E-01 

13500 

.116430E-01 

14000 

.126775E-01 

14500 

.125351E-01 

15000 

.127253E-01 

15500 

.129989E-01 

16000 

.128307E-01 

16500 

.128825E-01 

17000 

.128826E-01 

17500 

.126003E-01 

18000 

.126962E-01 

18500 

.125403E-01 

19000 

.124761E-01 

19500 

.122777E-01 

20000 

.128463E-01 

20500 

.128484E-01 

21000 

.127292E-01 

21500 

.125329E-01 

22000 

.123501E-01 

22500 

.123911E-01 

23000 

.123040E-01 

23500 

.121881E-01 

24000 

.122018E-01 

24500 

.121061E-01 

25000 

.121410E-01 

25500 

.120981E-01 

26000 

.120675E-01 

26500 

.121282E-01 

27000 

.119992E-01 

27500 

.119030E-01 

28000 

.118618E-01 

28500 

.117067E-01 

29000 

.117883E-01 

29500 

.118460E-01 

30000 

.117547E-01 

30500 

.118110E-01 
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31000  .117999E-01 
31500  .117988E-01 
32000  .117074E-01 
32500  .139721E-01 
33000  .142307E-01 
33500  .140911E-01 
34000  .139402E-01 
34500  .138819E-01 
35000  .137882E-01 
35500  .137614E-01 
36000  .136297E-01 
36500  .136593E-01 
37000  .136282E-01 
37500  .134956E~01 
38000  .134633E-01 
38500  .134175E-01 
39000  .133462E-01 
39500  .132598E-01 
40000  .131473E-01 
40500  .130197E-01 
41000  -129080E-01 
41500  .128247E-01 
42000  .128265E-01 
42500  .128098E--01 
43000  .127428E-01 
43500  .127139E-01 
44000  .126756E-01 
44500  .128702E-01 
45000  .128092E-01 
45500  .127654E-01 
46000  ■127176E-01 
46500  ■126982E-01 
47000  .126593E-01 
47500  .125463E-01 
48000  .124464E-01 
48500  .124563E-01 
49000  .123950E-01 
49500  .123531E-01 
50000  -123190E-01 


The  importance  failure  count  =  4035 
The  safety  factor  >  10  count  =  37 
The  reliability  index  =  .224703D+01 
The  binomial  statistics 

The  ave  of  Pu  (AvePu)  =  .123190D-01 
The  var  of  avePu  (var AvePu)  =  .243344D-06 
The  std  of  avePu  (stdAvePu)  =  .493299D-03 
The  cov  of  avePu  (covAvePu)  =  .4004390-01 


The  simulation  and  sample  mean  statistics 
The  ave  of  Pu  (avePu)  =  .1231900-01 

The  var  of  Pu  (varPu)b  =  .1376520+00 

The  var  of  Pu  (varPu)  =  .1376550+00 

The  std  of  Pu  (stdPu)  =  .3710190+00 

The  cov  of  Pu  (covPu)  =  .3011770+02 

The  var  of  avePu  (varAvePu)b  =  .2753050-05 

The  var  of  avePu  (varAvePu)  =  .2753100-05 

The  std  of  avePu  (stdAvePu)  =  .1659250-02 

The  cov  of  avePu  (covAvePu)  =  .1346910+00 

b  =  biased,  for  variance  computation  using  the  sample  size  N. 


Appendix  A  Case  Study  1 


A11 


Output  of  Importance  Sampling  -  IS02  {SFt  =1 .5) 

Title  =  Title  ASM-lbOl 

Total  number  of  random  variables  =  11 

*****  The  input  data  ***** 


X( 

X( 

X( 

X( 

X( 

X( 

X( 

X( 

X( 


rrv- 

— name- 

_ 

- mean - 

- Sigma - 

- cov - 

code 

Type 

1) 

STRU 

0 

1 

.1500000+00 

.2250000-01 

.1500000+00 

C 

NOR 

2) 

SOLT 

1 

3 

.2800000+02 

.8400000+01 

.3000000+00 

C 

NOR 

3) 

SOLT 

1 

5 

.1200000+00 

.1200000-01 

.1000000+00 

C 

NOR 

4) 

SORT 

1 

3 

.2800000+02 

.8400000+01 

.3000000+00 

C 

NOR 

5) 

SORT 

1 

5 

.1200000+00 

.1200000-01 

.1000000+00 

C 

NOR 

6) 

SOST 

0 

3 

.3000000+02 

.900000O+01 

,3000000+00 

C 

NOR 

7) 

WATR 

0 

1 

,5000000+01 

.1000000+01 

.2000000+00 

S 

NOR 

8) 

WATR 

0 

2 

.1500000+01 

.1000000+01 

.6666670+00 

S 

NOR 

9) 

VPLO 

1 

5 

.2000000+01 

.3000000+00 

.1500000+00 

C 

LOG 

10) 

VSLO 

1 

5 

.8000000+00 

.1200000+00 

.1500000+00 

C 

LOG 

11) 

VRLO 

1 

5 

.5500000+00 

.8250000-01 

.1500000+00 

C 

LOG 

The  limit  state  Z=:R-L=L( SF-1 )  or  Z=SF-1 


The 

The 

The 

The 

The 


simulation  cycle  =  50000 

initial  given  seed  =  -4175144 

sensitivity  test  ratio  = 

target  safety  factor  = 

initial  output  cycle  increment  = 


lOOOOOD+00 

150000D+01 

0 


«<  The  sensitivity  result  »> 

The  safety  factor  (SF)  -  ,182545E+01 


-x=rv-  — Name- 

— M=Mean — 

SF1:M*  .90 

SF2:M*1.10 

— SFl-SF - 

.-SF2-SF — 

Sense. — 

X (  1 )  STRU 

0 

1 

.15000+00 

.17960+01 

.18540+01 

-.28960-01 

.28860-01 

.90000+00 

X(  2)  SOLT 

1 

3 

.28000+02 

.17390+01 

.19150+01 

-.86800-01 

.89260-01 

.90000+00 

X(  3)  SOLT 

1 

5 

.12000+00 

.18780+01 

.17790+01 

.52980-01 

-.46150-01 

.11000+01 

X (  4 )  SORT 

1 

3 

.28000+02 

.17820+01 

.18720+01 

-.43900-01 

.46840-01 

.90000+00 

X(  5)  SORT 

1 

5 

.12000+00 

.17980+01 

.18530+01 

-.27580-01 

.27810-01 

.9000O+00 

X(  6)  SOST 

0 

3 

.30000+02 

.17430+01 

.19120+01 

-.82060-01 

.86330-01 

.90000+00 

X (  7 )  WATR 

0 

1 

.50000+01 

.18370+01 

.18140+01 

.11170-01 

-.11260-01 

.11000+01 

X (  8 )  WATR 

0 

2 

.15000+01 

.18310+01 

.18200+01 

.55010-02 

-.56180-02 

,11000+01 

X(  9)  VPLO 

1 

5 

.20000+01 

.18250+01 

.18250+01 

.OOOOO+OO 

.00000+00 

.10000+01 

X(  10)  VSLO 

1 

5 

.80000+00 

.17720+01 

.18800+01 

-.53660-01 

.54580-01 

.90000+00 

X(  11)  VRLO 

1 

5 

.55000+00 

.18760+01 

.17790+01 

.50200-01 

-.46110-01 

.llOOD+01 


«<  The  shifted  means  »> 


-x= 

=rv- 

— Name- 

— 

X( 

1) 

STRU 

0 

1 

X( 

2) 

SOLT 

1 

3 

X( 

3) 

SOLT 

1 

5 

X( 

4) 

SORT 

1 

3 

X( 

5) 

SORT 

1 

5 

X( 

6) 

SOST 

0 

3 

X( 

7) 

WATR 

0 

1 

X( 

8) 

WATR 

0 

2 

X( 

9) 

VPLO 

1 

5 

X( 

10) 

VSLO 

1 

5 

X( 

11) 

VRLO 

1 

5 

- M=Mean - 

.150000D+00 

.280000D+02 

.120000D+00 

.280000D+02 

.120000D+00 

.300000D+02 

.500000D+01 

.150000D+01 

.200000D+01 

.800000D+00 

.550000D+00 


Shifted. Mean 
.135000D+00 
.252000D+02 
.1320000+00 
.2520000+02 
.1080000+00 
.2700000+02 
.5500000+01 
.1650000+01 
.2000000+01 
.7200000+00 
.6050000+00 


- Ratio - 

.9000000+00 

.9000000+00 

.1100000+01 

.9000000+00 

.9000000+00 

.9000000+00 

.1100000+01 

.1100000+01 

.lOOOOOD+01 

.9000000+00 

.1100000+01 
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The  safety  factor  =  .143618D+01 

The  iteration  of  shifting  =  1 

*****  The  result  ***** 

The  revised  output  cycle  increment  =  500 

***  The  result  of  importance  sampling  simulation  *** 


No . of . cycle-  -ave . cum . Pu- 


500 

.114841E-01 

1000 

.110619E-01 

1500 

.115797E-01 

2000 

.940005E-02 

2500 

.882741E-02 

3000 

.820059E-02 

3500 

.888176E-02 

4000 

.832656E-02 

4500 

.826378E-02 

5000 

.889562E-02 

5500 

.968442E-02 

6000 

.933157E-02 

6500 

.934384E-02 

7000 

■918440E-02 

7500 

.903444E-02 

8000 

.916210E-02 

8500 

.911424E-02 

9000 

.940873E-02 

9500 

.958788E-02 

10000 

.102034E-01 

10500 

.102879E-01 

11000 

.997049E-02 

11500 

.101574E-01 

12000 

.102058E-01 

12500 

.996347E-02 

13000 

.976166E-02 

13500 

,988919E-02 

14000 

.101289E-01 

14500 

.100602E-01 

15000 

.100594E-01 

15500 

.994770E-02 

16000 

.982323E-02 

16500 

.966996E-02 

17000 

.961503E-02 

17500 

.973809E-02 

18000 

.952750E-02 

18500 

.948915E-02 

19000 

.958080E-02 

19500 

.953854E-02 

20000 

.961472E-02 

20500 

.955461E-02 

21000 

.958697E-02 

21500 

.945272E-02 

22000 

.935398E-02 

22500 

.924403E-02 

23000 

.912753E-02 

23500 

.935799E-02 

24000 

.934265E-02 

24500 

.933448E-02 

25000 

.931839E-02 

25500 

.942703E-02 

26000 

.950188E-02 

26500 

.943946E-02 

27000 

.937886E-02 

27500 

.925511E-02 

28000 

.919671E-02 

28500 

.921812E-02 

29000 

.910797E-02 

29500 

.924519E-02 

30000 

.925946E-02 

30500 

.932726E-02 
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31000 

31500 

32000 

32500 

33000 

33500 

34000 

34500 

35000 

35500 

36000 

36500 

37000 

37500 

38000 

38500 

39000 

39500 

40000 

40500 

41000 

41500 

42000 

42500 

43000 

43500 

44000 

44500 

45000 

45500 

46000 

46500 

47000 

47500 

48000 

48500 

49000 

49500 

50000 


.938460E-02 

.988154E-02 

.978469E-02 

.975847E-02 

.978069E-02 

.979664E-02 

.105081E-01 

.104348E-01 

.104052E-01 

.104833E-01 

.104445E-01 

.103646E-01 

.103946E-01 

.103583E-01 

.107490E-01 

.106795E-01 

.106835E-01 

.106057E‘-01 

.105313E-01 

.104907E-01 

.105045E-01 

.105086E-01 

.104417E-01 

.105012E-01 

.104821E-01 

.104270E-01 

.103719E-01 

.104500E-01 

.104161E-01 

.104109E-01 

.103741E-01 

.103418E-01 

.103337E-01 

.103191E-01 

.102739E-01 

.102667E-01 

■102907E-01 

.103235E-01 

.103543E-01 


The  importance  failure  count  =  4005 
The  safety  factor  >  10  count  =  31 
The  reliability  index  =  .231326D+01 
The  binomial  statistics 

The  ave  of  Pu  (AvePu)  =  .103543D-01 
The  var  of  avePu  (varAvePu)  =  .2049410-06 
The  std  of  avePu  (stdAvePu)  =  .4527040-03 
The  cov  of  avePu  (covAvePu)  =  .4372150-01 


The  simulation  and  sample  mean  statistics 
The  ave  of  Pu  (avePu)  =  .1035430-01 

The  var  of  Pu  (varPu)b  =  .2958340-01 

The  var  of  Pu  (varPu)  =  .2958400-01 

The  std  of  Pu  (stdPu)  =  .1720000+00 

The  cov  of  Pu  (covPu)  =  .1661150+02 

The  var  of  avePu  (varAvePu)b  =  .5916690-06 

The  var  of  avePu  (var AvePu)  -  .5916810-06 

The  std  of  avePu  (stdAvePu)  =  .7692080-03 

The  cov  of  avePu  (covAvePu)  =  .7428900—01 

b  =  biased,  for  variance  computation  using  the  sample  size  N. 
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Output  of  Direct  Monte  Carlo  Simulation  (MCS)  - 
OR01  (SFt  =1.82545) 

Title  =  Title  ASM-lbOl 

Total  number  of  random  variables  =  11 


*****  1 

-x=rv- 

The  input 

— name — 

data  ***** 

- mean - - 

- Sigma - 

- cov - 

code 

Type 

X( 

1) 

STRU 

0 

1 

.150000D+00 

.2250000-01 

.1500000+00 

C 

NOR 

X( 

2) 

SOLT 

1 

3 

.280000D+02 

.8400000+01 

.3000000+00 

C 

NOR 

X( 

3) 

SOLT 

1 

5 

.120000D+00 

.1200000-01 

.1000000+00 

C 

NOR 

X( 

4) 

SORT 

1 

3 

.280000D+02 

.8400000+01 

.3000000+00 

C 

NOR 

X( 

5) 

SORT 

1 

5 

.120000D+00 

.1200000-01 

.1000000+00 

C 

NOR 

X( 

6) 

SOST 

0 

3 

.300000D+02 

.9000000+01 

.3000000+00 

C 

NOR 

X( 

7) 

WATR 

0 

1 

.500000D+01 

.1000000+01 

.2000000+00 

S 

NOR 

X( 

8) 

WATR 

0 

2 

.1500000+01 

.1000000+01 

.6666670+00 

S 

NOR 

X( 

9) 

VPLO 

1 

5 

.2000000+01 

.3000000+00 

.1500000+00 

C 

LOG 

X( 

10) 

VSLO 

1 

5 

.800000D+00 

.1200000+00 

.1500000+00 

c 

LOG 

X( 

11) 

VRLO 

1 

5 

.5500000+00 

.8250000-01 

.1500000+00 

c 

LOG 

The 

1  limit  state 

2=R-L=L(SF-1) 

or  Z=SF-1 

The  simulation  cycle  =  50000 

The  initial  given  seed  =  -362646 

The  sensitivity  test  ratio  = 

The  target  safety  factor  = 

The  initial  output  cycle  increment  = 


lOOOOOD+00 

lOOOOOD+00 

0 


*****  The  result  ***** 

The  revised  output  cycle  increment  =  500 

***  The  result  of  direct  simulation  *** 


No. of. cycle-  -ave.cum.Pu- 


500 

.400000E-02 

1000 

.130000E-01 

1500 

.113333E-01 

2000 

.lOOOOOE-01 

2500 

.108000E-01 

3000 

.103333E-01 

3500 

.111429E-01 

4000 

.107500E-01 

4500 

.104444E-01 

5000 

.102000E-01 

5500 

.lOOOOOE-01 

6000 

.lOOOOOE-01 

6500 

.953846E-02 

7000 

.101429E-01 

7500 

.986667E-02 

8000 

.lOOOOOE-01 

8500 

.101176E-01 

9000 

.lOllllE-01 

9500 

.989474E-02 

10000 

.980000E-02 

10500 

.100952E-01 

11000 

.lOOOOOE-01 

11500 

.100870E-01 

12000 

.101667E-01 

12500 

.100800E-01 

13000 

.992308E-02 

13500 

.985185E-02 

14000 

.992857E-02 

14500 

.lOOOOOE-01 

15000 

.100667E-01 

15500 

.100645E-01 

16000 

.101250E-01 

16500 

.lOOOOOE-01 
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17000 

.100588E-01 

17500 

.994286E-02 

18000 

.988889E-02 

18500 

.101622E-01 

19000 

.100526E~01 

19500 

.994872E-02 

20000 

.995000E-02 

20500 

.995122E-02 

21000 

.995238E-02 

21500 

.lOOOOOE-01 

22000 

.995455E-02 

22500 

.102667E-01 

23000 

.101739E-01 

23500 

.101702E-01 

24000 

.102083E-01 

24500 

.102041E-01 

25000 

.102800E~01 

25500 

.101176E-01 

26000 

.101154E-01 

26500 

.101509E-01 

27000 

.101111E~01 

27500 

.101091E-01 

28000 

.996429E~02 

28500 

.989474E-02 

29000 

.982759E-02 

29500 

.983051E-02 

30000 

.986667E-02 

30500 

.996721E~02 

31000 

.987097E-02 

31500 

.971429E-02 

32000 

.971875E-02 

32500 

.969231E-02 

33000 

.963636E-02 

33500 

.958209E-02 

34000 

.955882E-02 

34500 

.959420E-02 

35000 

.948571E-02 

35500 

.949296E-02 

36000 

.955556E-02 

36500 

.950685E-02 

37000 

.951351E-02 

37500 

.952000E-02 

38000 

.963158E-02 

38500 

.961039E-02 

39000 

.958974E-02 

39500 

.954430E-02 

40000 

.955000E-02 

40500 

.958025E-02 

41000 

.953659E-02 

41500 

.949398E-02 

42000 

.961905E-02 

42500 

.969412E-02 

43000 

.965116E-02 

43500 

.967816E-02 

44000 

.965909E-02 

44500 

.966292E-02 

45000 

.966667E-02 

45500 

.971429E-02 

46000 

.967391E-02 

46500 

.965591E~02 

47000 

.965957E-02 

47500 

.972632E-02 

48000 

.970833E-02 

48500 

.967010E-02 

49000 

.963265E-02 

49500 

.959596E-02 

50000 

.960000E-02 

The  importance  failure  count  = 
The  safety  factor  >  10  count  = 


480 
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The  reliability  index 

The  binomial  statistics 
The  ave  of  Pu  (AvePu) 
The  var  of  avePu  (varAvePu) 
The  std  of  avePu  (stdAvePu) 
The  cov  of  avePu  (covAvePu) 


.234162D+01 


.960000D-02 

.190157D-06 

.436070D-03 

.454239D-01 
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Output  of  Direct  Monte  Carlo  Simulation  (MCS)  - 
DR02  (SF,  =1.82545) 

Title  =  Title  ASM-lbOl 

Total  number  of  random  variables  =  11 

■k****  The  input  data  ***** 


-x= 

=rv- 

X( 

1) 

X{ 

2) 

X{ 

3) 

X( 

4) 

X( 

5) 

X( 

6) 

X( 

7) 

X( 

8) 

X( 

9) 

X( 

10) 

X( 

11) 

— name- 

_ 

- mean - 

- Sigma - 

- cov - 

code 

Type 

STRU 

0 

1 

.150000D+00 

.225000D-01 

.1500000+00 

C 

NOR 

SOLT 

1 

3 

.280000D+02 

.8400000+01 

.3000000+00 

C 

NOR 

SOLT 

1 

5 

.120000D+00 

.120000D-01 

.1000000+00 

C 

NOR 

SORT 

1 

3 

.280000D+02 

.8400000+01 

.3000000+00 

C 

NOR 

SORT 

1 

5 

.120000D+00 

.1200000-01 

.1000000+00 

C 

NOR 

SOST 

0 

3 

.300000D+02 

.900000O+01 

.3000000+00 

C 

NOR 

WATR 

0 

1 

.500000D+01 

.1000000+01 

,2000000+00 

S 

NOR 

WATR 

0 

2 

.150000D+01 

.1000000+01 

.6666670+00 

S 

NOR 

VPLO 

1 

5 

.200000D+01 

.3000000+00 

,1500000+00 

C 

LOG 

VSLO 

1 

5 

.800000D+00 

.1200000+00 

.1500000+00 

C 

LOG 

VRLO 

1 

5 

.550000D+00 

.8250000-01 

.1500000+00 

C 

LOG 

The  limit  state  Z=R-L=L( SF-1 )  or  Z=SF-1 


The  simulation  cycle  =  50000 

The  initial  given  seed  =  -4175144 

The  sensitivity  test  ratio  = 

The  target  safety  factor  = 

The  initial  output  cycle  increment  = 


lOOOOOD+00 

lOOOOOD+00 

0 


*****  The  result  ***** 

The  revised  output  cycle  increment  -  500 

***  The  result  of  direct  simulation  *** 


No. of. cycle-  -ave.cum.Pu- 


500  .800000E-02 

1000  .130000E-01 

1500  -lOOOOOE-01 
2000  .950000E-02 

2500  .104000E-01 

3000  .966667E-02 

3500  .971429E-02 

4000  .112500E-01 

4500  .108889E-01 

5000  .108000E-01 

5500  .107273E-01 

6000  .llOOOOE-01 
6500  .107692E-01 

7000  .111429E-01 

7500  .113333E-01 

8000  .111250E-01 

8500  .107059E-01 

9000  .108889E-01 

9500  .105263E-01 

10000  .108000E-01 

10500  .109524E-01 

11000  .106364E-01 

11500  .106087E-01 

12000  .105833E-01 

12500  .103200E-01 

13000  .103077E-01 

13500  .102963E-01 

14000  .102857E-01 

14500  .103448E-01 

15000  .104667E-01 

15500  .105161E-01 

16000  .104375E-01 

16500  .106061E-01 
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17000 

.105882E~01 

17500 

.107429E-01 

18000 

.106111E-01 

18500 

.105946E-01 

19000 

.107895E-01 

19500 

.107692E-01 

20000 

.107500E-01 

20500 

.107317E-01 

21000 

.105714E-01 

21500 

.105581E-01 

22000 

.107273E-01 

22500 

.107111E-01 

23000 

.106087E-01 

23500 

.105957E-01 

24000 

.105833E-01 

24500 

.106122E-01 

25000 

.106800E-01 

25500 

.106275E-01 

26000 

.106923E-01 

26500 

.105660E-01 

27000 

.106667E-01 

27500 

.105455E-01 

28000 

.105357E-01 

28500 

.104912E-01 

29000 

.104483E-01 

29500 

.103729E-01 

30000 

.103667E-01 

30500 

.104590E-01 

31000 

.104516E-01 

31500 

.103492E-01 

32000 

.102500E-01 

32500 

.102154E-01 

33000 

.101818E-01 

33500 

.100597E-01 

34000 

.lOOOOOE-01 

34500 

.100580E-01 

35000 

.994286E-02 

35500 

.991549E-02 

36000 

.986111E-02 

36500 

.980822E-02 

37000 

.978378E-02 

37500 

.973333E-02 

38000 

.971053E-02 

38500 

.971429E-02 

39000 

.966667E-02 

39500 

.969620E-02 

40000 

.975000E-02 

40500 

.975309E-02 

41000 

.970732E-02 

41500 

.973494E-02 

42000 

.976190E~02 

42500 

.981176E-02 

43000 

.981395E-02 

43500 

.981609E-02 

44000 

.977273E~02 

44500 

.979775E-02 

45000 

.100222E-01 

45500 

.100440E-01 

46000 

.997826E-02 

46500 

.993548E-02 

47000 

.991489E-02 

47500 

.989474E-02 

48000 

.987500E-02 

48500 

.987629E-02 

49000 

.989796E-02 

49500 

.989899E-02 

50000 

.lOOOOOE-01 

The  importance  failure  count  = 
The  safety  factor  >  10  count  = 


500 
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The  reliability  index 


The 

The 

The 

The 

The 


binomial  statistics 
ave  of  Pu  (AvePu) 
var  of  avePu  (var AvePu) 
std  of  avePu  (stdAvePu) 
cov  of  avePu  (covAvePu) 


.232635D+01 


.lOOOOOD-01 

.1980000-06 

.4449720-03 

.4449720-01 
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Appendix  B 

Input  and  Output  Data  for  Case 
Study  2  (Correlation) 


Input  Data  for  CSLIDE 

001  TITL  RETAINING  WALL  with  correlation 
002  TITL 
003  TITL 
004  TITL 

005  STRU  8  1.500000E-01  O.OOOOOOE+00  l.OOOOOOE+00 

006  0.00  0.00 

007  0.00  2.00 

008  6.00  2.00 

009  6.00  14.00 

010  8.00  14.00 

011  8.00  2.00 

012  12.00  2.00 

013  12.00  0,00 

014  SOLT  1  1  .280000E+02  .500000E-01  .120000E+00  .140000E+02 

015  -500.00  14.00 

016  SORT  1  1  .280000E+02  .500000E-01  .120000E+00  .400000E+01 

017  500.00  4.00 

018  SOST  3.000000E+01  0.500000E-01 

019  METH  1 

020  WATR  5.000000E+00  1.500000E+00  6.250000E-02  -1 

021  END 


Input  Data  for  Reliability  Methods 

TITL  ASM/MCS  data  file  for  X0075dle.dat  -  correlation 

STRUOOOl  1.500000E-01  1.500000E-01  C  NOR 

SOLT0103  2.800000E+01  3.G00000E-01  C  NOR 

SOLT0104  5.000000E-02  2.000000E-01  C  LOG  -0.5 

SOLT0105  1.200000E-01  l.OOOOOOE-01  C  NOR 

SORT0103  2.800000E+01  3.000000E-01  C  NOR 

SORT0104  5.000000E-02  2.000000E-01  C  LOG  -0.5 

SORT0105  1.200000E-01  l.OOOOOOE-01  C  NOR 

SOST0003  3.000000E+01  3.000000E-01  C  NOR 

SOST0004  5.000000E-02  l.OOOOOOE-02  C  NOR  -0.5 

WATROOOl  5.000000E+00  l.OOOOOOE+00  S  NOR 

WATR0002  l.OOOOOOE+00  l.OOOOOOE+OO  S  NOR 

END 
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Output  of  Advanced  Second  Moment  (ASM) 
Method 


Title  =  ASM/MCS  data  file  for  X0075dle.dat  -  correlation 
Total  number  of  random  variables  =  11 


*****  The  input  data  ***** 


-x= 

=rv- 

— name- 

_ 

- mean - 

- sigma - 

- cov - 

code 

Type 

X( 

1) 

STRU 

0 

1 

.1500000+00 

.2250000-01 

.1500000+00 

C 

NOR 

X( 

2) 

SOLT 

1 

3 

.2800000+02 

.8400000+01 

.300000O+00 

C 

NOR 

X( 

3) 

SOLT 

1 

4 

.5000000-01 

.1000000-01 

-2000000+00 

C 

NOR 

Correlation 

for 

SOLT  1  =  -. 

500000O+00 

X( 

4) 

SOLT 

1 

5 

.1200000+00 

.1200000-01 

.1000000+00 

C 

NOR 

X 

5) 

SORT 

1 

3 

.2800000+02 

.8400000+01 

.3000000+00 

C 

NOR 

X 

6) 

SORT 

1 

4 

.5000000-01 

.1000000-01 

.2000000+00 

C 

NOR 

Correlation 

for 

SORT  1  =  -. 

5000000+00 

X( 

7) 

SORT 

1 

5 

.1200000+00 

.1200000-01 

.1000000+00 

C 

NOR 

X{ 

8) 

SOST 

0 

3 

.3000000+02 

.9000000+01 

.3000000+00 

C 

NOR 

X 

9) 

SOST 

0 

4 

.5000000-01 

.5000000-03 

.1000000-01 

C 

NOR 

Correlation 

for 

SOST  0  =  -. 

5000000+00 

X( 

10) 

WATR 

0 

1 

.5000000+01 

.1000000+01 

.2000000+00 

S 

NOR 

X( 

11) 

WATR 

0 

2 

.1000000+01 

.1000000+01 

.1000000+01 

S 

NOR 

The  tolerance  of  reliability  index  (beta)  =  .lOOOOOD-01 

The  ratio  of  dx/x  for  random  variables  =  .lOOOOOD-02 

The  max.  iterations  for  finding  beta  =  10 

The  increment  for  finding  beta  =  .2500 

The  l=absolute  2=relative  tolerance  for  beta  =  2 


The  limit  state  Z=R-L=L( SF-1 )  or  Z=SF-1 

*****  The  result  ***** 

The  iteration  No.  =1 
-i - meanEN - - stdEN - 

1  .150000D+00  .225000D-01 

2  .280000D+02  .840000D+01 

3  .500000D-01  -lOOOOOD-01 

4  .120000D+00  .120000D~01 

5  .280000D+02  .840000D+01 

6  .500000D-01  .lOOOOOD-01 

7  .120000D+00  .120000D-01 

8  .300000D+02  .9000000+01 

9  .500000D-01  .5000000-03 

10  .500000O+01  .lOOOOOO+Ol 

11  .1000000+01  .1000000+01 

The  reliability  index  (beta) 

The  iteration  No.  =  2 
-i - meanEN - - stdEN - 

1  .1500000+00  .2250000-01 

2  .2800000+02  .8400000+01 

3  .5000000-01  .1000000-01 
4  .1200000+00  .1200000-01 

5  -2800000+02  .8400000+01 

6  .5000000-01  .1000000-01 
7  .1200000+00  .1200000-01 

8  .3000000+02  .9000000+01 

9  .5000000-01  .5000000-03 

10  .5000000+01  .1000000+01 
11  .1000000+01  .1000000+01 

The  reliability  index  (beta)  =  .2767850+01 

The  unsatisfactory  performance  probability  =  .2821360-02 


- alpha - - dp - -Partial.  SF- 

.1569510+00  .1402630+00  .9350890+00 

.3092180+00  .1554510+02  .5551820+00 

.4424360+00  .5630160-01  .1126030+01 

-.2065910+00  .1268350+00  .1056960+01 

.8814810-01  .2482990+02  .8867830+00 

.1071580+00  .5134350-01  .1026870+01 

.9441390-01  .1168760+00  .9739680+00 

.3860140+00  .1089860+02  -3632860+00 

.6659830+00  .5052900-01  .1010580+01 

-.8176050-01  .5225430+01  .1045090+01 

-.1180170+00  .1325390+01  .1325390+01 

=  .2757180+01 


- alpha - - dp - -Partial.  SF- 

.7326000-01  .1454380+00  .9695840+00 

.2981470+00  .1505480+02  .5376730+00 

.4707810+00  .5715860-01  .1143170+01 

-.2028870+00  .1267390+00  .1056160+01 

.1676320+00  .2020510+02  .7216120+00 

.2903460+00  .5463980-01  . 1092800+01 

.7686910-01  . 1174470+00  .9787240+00 

.3620250+00  .1196340+02  .3987810+00 

.6270460+00  .5050100-01  .1010020+01 

-.1960540-01  .5054260+01  .1010850+01 

. 0000000+00  . 1000000+01  .1000000+01 
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Output  of  Advanced  Second  Moment  (ASM) 
Method  (Without  Correlation) 

Title  =  ASM/MCS  data  file  for  X0075dld.dat  -  Without  correlation 
Total  number  of  random  variables  =  11 


*****  The  input  data  ***** 


-x= 

XV- 

— name- 

— .  . 

- mean - 

- Sigma - 

- cov - 

code 

Type 

X( 

1) 

STRU 

0 

1 

.1500000+00 

.2250000-01 

.1500000+00 

c 

NOR 

X( 

2) 

SOLT 

1 

3 

.2800000+02 

.8400000+01 

.3000000+00 

c 

NOR 

X( 

3) 

SOLT 

1 

4 

.5000000-01 

.1000000-01 

.2000000+00 

c 

NOR 

X( 

4) 

SOLT 

1 

5 

.1200000+00 

.1200000-01 

.1000000+00 

c 

NOR 

X( 

5) 

SORT 

1 

3 

.2800000+02 

.8400000+01 

.3000000+00 

c 

NOR 

X( 

6) 

SORT 

1 

4 

.5000000-01 

.1000000-01 

.2000000+00 

c 

NOR 

X( 

7) 

SORT 

1 

5 

.1200000+00 

.1200000-01 

.1000000+00 

c 

NOR 

X( 

8) 

SOST 

0 

3 

.3000000+02 

.9000000+01 

.3000000+00 

c 

NOR 

X( 

9) 

SOST 

0 

4 

.5000000-01 

.5000000-03 

.1000000-01 

c 

NOR 

X( 

10) 

WATR 

0 

1 

.5000000+01 

.1000000+01 

.2000000+00 

s 

NOR 

X( 

11) 

WATR 

0 

2 

.1000000+01 

.1000000+01 

.1000000+01 

s 

NOR 

The 

tolerance 

of 

reliability 

index  ( beta ) 

=  .1000000 

-01 

The  ratio  of  dx/x  for  random  variables  =  .lOOOOOD-02 

The  max.  iterations  for  finding  beta  =  10 

The  increment  for  finding  beta  =  .2500 

The  l=absolute  2=relative  tolerance  for  beta  =  2 

The  limit  state  Z=R-L:=L(SF~1)  or  Z^SF-1 

*****  The  result  ***** 


The  iteration  No.  .=  1 

-i - meanEN - - stdEN - - alpha - - dp - -Partial.  SF- 


1  .150000D+00  .225000D-01  .154239D+00  .140666D+00  .937776D+00 

2  .280000D+02  .840000D+01  .554903D+00  .154637D+02  .552275D+00 

3  .500000D-01  .lOOOOOD-01  .528479D-01  .485787D-01  .971573D+00 

4  .120000D+00  .120000D-01  -.203022D+00  .126552D+00  .1054600+01 

5  .2800000+02  .8400000+01  .1474240+00  .2466940+02  .8810500+00 

6  .5000000-01  .1000000-01  .2582610-01  .4930540-01  .9861080+00 

7  .1200000+00  .1200000-01  .9278280-01  .1170060+00  .9750460+00 

8  .3000000+02  .9000000+01  .7572080+00  .1167130+02  .3890450+00 

9  .5000000-01  .5000000-03  .1482610-02  .4999800-01  .9999600+00 

10  .5000000+01  .1000000+01  -.8034800-01  .5216100+01  .1043220+01 

11  .1000000+01  .1000000+01  -.1159780+00  .1311930+01  .1311930+01 

The  reliability  index  (beta)  =  .2689510+01 

The  iteration  No.  -2 

-i - meanEN - - stdEN - - alpha - - dp - -Partial -SF- 

1  .1500000+00  .2250000-01  .5010940-01  .1469780+00  .9798560+00 

2  .2800000+02  .8400000+01  .6444040+00  .1349310+02  .4818950+00 

3  .5000000-01  .1000000-01  .3645870-01  .4902290-01  .9804580+00 

4  .1200000+00  .1200000-01  -.1674580+00  .1253850+00  .1044880+01 

5  .2800000+02  .8400000+01  .2834560+00  .2161880+02  .7721000+00 

6  .5000000-01  .1000000-01  .2793440-01  .4925140-01  .9850270+00 

7  .1200000+00  .1200000-01  .9240950-01  .1170280+00  .9752340+00 

8  .3000000+02  .9000000+01  .6744100+00  .1373310+02  .4577700+00 

9  .5000000-01  .5000000-03  .1363660-02  .4999820-01  .9999630+00 

10  .5000000+01  .1000000+01  -.2888410-01  .5077410+01  .1015480+01 

11  .1000000+01  .1000000+01  -.8676510-01  .1232530+01  -1232530+01 

The  reliability  index  (beta)  =  .2680020+01 

The  unsatisfactory  performance  probability  =  .3680900-02 
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Output  of  Direct  Monte  Carlo  Simulation 

Title  =  ASM/MCS  data  file  for  X0075dle.dat  -  correlation 
Total  number  of  random  variables  =  11 


*****  The  input  data  ***** 


-x= 

=rv- 

— name- 

_ 

- mean - 

- Sigma - 

- cov - 

code 

Type 

X( 

1) 

STRU 

0 

1 

.150000D+00 

.225000D-01 

.150000O+00 

C 

NOR 

X 

2) 

SOLT 

1 

3 

.280000D+02 

.840000D+01 

.3000000+00 

C 

NOR 

X 

3) 

SOLT 

1 

4 

.500000D-01 

.lOOOOOD-01 

.2000000+00 

C 

NOR 

Correlation 

for 

SOLT  1  = 

500000D+00 

X( 

4) 

SOLT 

1 

5 

.120000D+00 

.120000D-01 

-1000000+00 

C 

NOR 

X( 

5) 

SORT 

1 

3 

.280000D+02 

.840000D+01 

.3000000+00 

C 

NOR 

X( 

6) 

SORT 

1 

4 

.500000D-01 

.lOOOOOD-01 

.200000O+00 

C 

NOR 

Correlation 

for 

SORT  1  =  -. 

500000D+00 

X( 

7) 

SORT 

1 

5 

.120000D+00 

.120000D-01 

.1000000+00 

C 

NOR 

X( 

8) 

SOST 

0 

3 

.300000D+02 

.9000000+01 

.300000O+00 

C 

NOR 

X 

9) 

SOST 

0 

4 

.500000D-01 

.5000000-03 

.1000000-01 

C 

NOR 

Correlation 

for 

SOST  0  =  “. 

5000000+00 

X( 

10) 

WATR 

0 

1 

.500000D+01 

-1000000+01 

.2000000+00 

S 

NOR 

x| 

11) 

WATR 

0 

2 

.lOOOOOD+01 

.1000000+01 

.1000000+01 

S 

NOR 

The  limit  state  Z=R-L=L( SF-1 )  or  Z=SF-1 
The  simulation  cycle  =  3000 

The  initial  given  seed  =  0 

The  simulation  output  increment  =  10 

The  target  safety  factor  =  .180000D+01 

The  safety  factor  shift  ratio  -  .lOOOOOD+00 

The  safety  factor  (SF)  -  .193198E+01 

*****  The  result  ***** 

The  initial  random  seed  =  -1142561 

***  The  result  of  direct  simulation  *** 


No. of. cycle-  -ave.cum.Pf-  -COV(avePf)-  B.COV(avePf) 


10  .OOOOOOE+00 
20  .OOOOOOE+00 

30  .333333E-01  .OOOOOOE+00  .983192E+00 

40  .250000E-01  .OOOOOOE+00  .987421E+00 

50  .200000E-01  .OOOOOOE+00  .989949E+00 

60  .166667E-01  .OOOOOOE+00  .991632E+00 

70  . 142857E-01  .OOOOOOE+00  .992831E+00 

80  .125000E-01  .OOOOOOE+00  .993730E+00 

90  . llllllE-01  .OOOOOOE+00  .994429E+00 

100  . lOOOOOE-01  .OOOOOOE+00  .994987E+00 

500  .600000E-02  .OOOOOOE+00  .575616E+00 

1000  .700000E-02  .OOOOOOE+00  .376639E+00 

1500  . 733333E-02  . OOOOOOE+00  . 300404E+00 

2000  .600000E-02  .OOOOOOE+00  .287808E+00 

2500  .560000E-02  . OOOOOOE+00  .266512E+00 

2900  .517241E-02  .OOOOOOE+00  .257530E+00 

2910  .515464E-02  .OOOOOOE+00  .257533E+00 

2920  .513699E-02  .OOOOOOE+00  .257535E+00 

2930  .511945E-02  .OOOOOOE+00  .257537E+00 

2940  .510204E-02  .OOOOOOE+00  .257539E+00 

2950  .508475E-02  .OOOOOOE+00  .257542E+00 

2960  .506757E-02  .OOOOOOE+00  .257544E+00 

2970  .505051E-02  .OOOOOOE+00  .257546E+00 

2980  .503356E-02  .OOOOOOE+00  .257548E+00 

2990  .501672E-02  .OOOOOOE+00  .257550E+00 

3000  .500000E-02  .OOOOOOE+00  .257553E+00 


The  importance  failure  count  =  15 
The  safety  factor  >10  count  =  0 
The  reliability  index  (beta)  =  .257583E+01 
The  binomial  statistics 

The  ave  of  Pf  (AvePf)  =  .500000D-02 
The  var  of  avePf  (varAvePf)  =  .165833D-05 
The  std  of  avePf  (stdAvePf)  =  .128776D-02 
The  cov  of  avePf  (covAvePf)  =  . 25755 3D+ 00 
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Output  of  Importance  Sampling  (SFt  =1 .8) 

Title  =  ASM/MCS  data  file  for  X0075dle.dat  -  correlation 


Total  number  of 

random  variables  =  11 

*****  The  input 

data  ***** 

-x= 

rv-  — name — 

- mean - - 

- Sigma - 

- cov - 

code 

Type 

X( 

1)  STRU 

0  1 

.150000D+00 

.2250000-01 

.1500000+00 

C 

NOR 

X( 

2)  SOLT 

1  3 

,2800000+02 

.8400000+01 

.3000000+00 

c 

NOR 

X( 

3)  SOLT 

1  4 

.5000000-01 

.1000000-01 

.2000000+00 

c 

NOR 

Correlation 

for 

SOLT  1  =  -.5000000+00 

X( 

4)  SOLT 

1  5 

.1200000+00 

.1200000-01 

.1000000+00 

c 

NOR 

X{ 

5)  SORT 

1  3 

.2800000+02 

.8400000+01 

.3000000+00 

c 

NOR 

X( 

6)  SORT 

1  4 

.5000000-01 

.1000000-01 

.2000000+00 

c 

NOR 

Correlation 

for 

SORT  1  =  -.5000000+00 

X( 

7)  SORT 

1  5 

.1200000+00 

.1200000-01 

-1000000+00 

c 

NOR 

X( 

8)  SOST 

0  3 

.3000000+02 

.9000000+01 

.3000000+00 

c 

NOR 

X( 

9)  SOST 

0  4 

.5000000-01 

.5000000-03 

.1000000-01 

c 

NOR 

Correlation 

for 

SOST  0  =  -.5000000+00 

X( 

10)  WATR 

0  1 

.5000000+01 

.1000000+01 

.2000000+00 

s 

NOR 

X( 

11)  WATR 

0  2 

.1000000+01 

.1000000+01 

.1000000+01 

s 

NOR 

The 

limit  state 

Z=R-L=L(SF-1) 

or  Z=SF-1 

The 

simulation  cycle  = 

3000 

The 

initial 

given  seed  = 

0 

The 

simulation  output  increment  = 

10 

The 

target  safety  factor 

^  .1800000+01 

The 

safety  factor  shift  ratio 

=  .1000000+00 

«<  The  sensitivity  result  »> 

The  safety  factor  (SF)  -  .193198E+01 


-x=rv-  — Name- 

Sense. — 

— 

— M=Mean — 

SF1:M*  .90 

SF2:M*1.10 

— SFl-SF-  ■ 

— SF2-SF— 

X(  1)  STRU 
.90000+00 

0 

1 

.15000+00 

.18890+01 

.19750+01 

-.43120-01 

.42930-01 

X(  2)  SOLT 
.90000+00 

1 

3 

.28000+02 

.18550+01 

.20100+01 

-.76510-01 

.78420-01 

X(  3)  SOLT 
.90000+00 

1 

4 

.50000-01 

.19210+01 

.19430+01 

-.11030-01 

.11050-01 

X(  4)  SOLT 
.11000+01 

1 

5 

.12000+00 

.20270+01 

.18550+01 

.94750-01 

-.76920-01 

X(  5)  SORT 
.90000+00 

1 

3 

.28000+02 

.19120+01 

.19540+01 

-.19690-01 

.21560-01 

X(  6)  SORT 
.90000+00 

1 

4 

.50000-01 

.19270+01 

.19370+01 

-.54370-02 

.54270-02 

X(  7)  SORT 
.90000+00 

1 

5 

.12000+00 

.18930+01 

.19710+01 

-.38560-01 

.39030-01 

X(  8)  SOST 
-9000O+00 

0 

3 

.30000+02 

.18290+01 

.20400+01 

-.10310+00 

.10840+00 

X(  9)  SOST 
.90000+00 

0 

4 

.50000-01 

.19260+01 

.19380+01 

-.62110-02 

.62080-02 

X(  10)  WATR 
.11000+01 

0 

1 

.50000+01 

.19490+01 

.19150+01 

.16700-01 

-.16840-01 

X(  11)  WATR 

0 

2 

.10000+01 

.19620+01 

.19520+01 

.29780-01 

.20080-01 

.llOOD+01 

«<  The  shifted  means  »> 

**NOTE**  correlated  random  variables  were  not  shifted** 


-x= 

=rv- 

— Name- 

- M=Mean - 

Shifted. Mean 

- Ratio - 

X( 

1) 

STRU 

0 

1 

.1500000+00 

.1350000+00 

.9000000+00 

X( 

2) 

SOLT 

1 

3 

.2800000+02 

.2800000+02 

.1000000+01 

X( 

3) 

SOLT 

1 

4 

.5000000-01 

.5000000-01 

.1000000+01 

X( 

4) 

SOLT 

1 

5 

.1200000+00 

.1320000+00 

.1100000+01 

X( 

5) 

SORT 

1 

3 

.2800000+02 

.2800000+02 

.1000000+01 

X( 

6) 

SORT 

1 

4 

.5000000-01 

.5000000-01 

.1000000+01 

X( 

7) 

SORT 

1 

5 

.1200000+00 

.1080000+00 

.9000000+00 

X( 

8) 

SOST 

0 

3 

.3000000+02 

.3000000+02 

.1000000+01 

X( 

9) 

SOST 

0 

4 

.5000000-01 

.5000000-01 

.1000000+01 

X( 

10) 

WATR 

0 

1 

.5000000+01 

.5500000+01 

.1100000+01 

X( 

11) 

WATR 

0 

2 

.1000000+01 

.1100000+01 

.1100000+01 

The  safety  factor  =  .158121D+01 

The  iteration  of  shifting  ==  1 
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*****  The  result  ***** 

The  initial  random  seed  =  -9527654 

***  The  result  of  importance  sampling  simulation  *** 

No. of. cycle-  -ave.cum.Pf-  -COV(avePf )-  B-COV(avePf) 


10  .OOOOOOE+00 
20  .OOOOOOE+00 
30  .OOOOOOE+00 
40  .300743E-03 

50  .240595E-03 

60  .200495E-03 

70  .171853E-03 

80  .150372E-03 

90  .133664E-03 

100  .120297E-03 

500  .333459E-02 

1000  .199691E-02 

1500  .207756E-02 

2000  .355221E-02 

2500  .296705E-02 

2900  .292495E-02 

2910  .293193E-02 

2920  .292189E-02 

2930  .291191E-02 

2940  .290201E-02 

2950  .289217E-02 

2960  .288240E-02 

2970  .289210E-02 

2980  .288240E-02 

2990  .287276E-02 

3000  .286318E-02 


.698519E+00 

.700387E+00 

.701619E+00 

.702492E+00 

.703144E+00 

.703649E+00 

.704051E+00 

.991259E+00 

.835725E+00 

.594589E+00 

.495030E+00 

.475451E+00 

.423550E+00 

.421126E+00 

.421127E+00 

.421128E+00 

.421129E+00 

.421130E+00 

.421131E+00 

.418355E+00 

.418356E+00 

.418357E+00 

.418358E+00 


.911605E+01 

.911633E+01 

.911651E+01 

.911664E+01 

.911674E+01 

.911681E+01 

.911687E+01 

.773159E+00 

.706946E+00 

.565882E+00 

.374509E+00 

.366625E+00 

.342851E+00 

.341853E+00 

.341855E+00 

.341856E+00 

.341858E+00 

.341860E+00 

.341861E+00 

.340711E+00 

.340712E+00 

.340714E+00 

.340716E+00 


The  importance  failure  count  =  33 
The  safety  factor  >  10  count  =  0 
The  reliability  index  (beta)  =  .276305E+01 


The  binomial  statistics 


The 

ave 

of 

Pf 

(AvePf) 

= 

.286318D-02 

The 

var 

of 

avePf 

( var AvePf ) 

= 

.951662D-06 

The 

std 

of 

avePf 

( stdAvePf ) 

= 

.975532D-03 

The 

cov 

of 

avePf 

(covAvePf ) 

-340716D+00 

The 

simulation  and  sample  mean 

statistics 

The 

ave 

of 

Pf 

(avePf ) 

.2863180-02 

The 

var 

of 

Pf 

(varPf )b 

= 

.4303000-02 

The 

var 

of 

Pf 

(varPf ) 

= 

.4304430-02 

The 

std 

of 

Pf 

(stdPf ) 

= 

.6560820-01 

The 

cov 

of 

Pf 

(covPf ) 

= 

.2291440+02 

The 

var 

of 

avePf 

(varAvePf )b 

.1434330-05 

The 

var 

of 

avePf 

( var AvePf ) 

= 

.1434810-05 

The 

std 

of 

avePf 

( StdAvePf ) 

= 

.1197840-02 

The 

cov 

of 

avePf 

(covAvePf ) 

= 

.4183580+00 

b  =  biased,  for  variance  computation  using  the  sample  size  N 
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Appendix  C 
Notation 


a 

ASM 

b 

bo 

bi 

c 

CDF 

COV 

Cov 

CSLIDE 

d 

DIFF 

DR 

dx 

f 

F 

FOSM 

G 

h 

I 

INV 

IS 

L 

LRFD 

MCS 

MOM 

N 

P 

PAR 

PDF 

R 

RCSLIDE 

SF 


Value  of  a  variable  X 
Advanced  second  moment  (method) 

Value  of  a  variable  X  or  regression  coefficient 

Intercept  of  a  regression  line  between  Zi  dXidXi 

Slope  of  the  regression  line 

Value  of  a  variable  X 

Cumulative  distribution  function 

Coefficient  of  variation 

Covariance 

Software  for  stability  assessment  of  concrete  gravity  structures 

Total  derivative 

Difference 

Direct  simulation 

Differential  of  x 

Probability  density  function 

Cumulative  distribution  function  or  force 

First-order  second  moment  method 

Performance  function 

Importance  density  function 

Performance  indicator 

Inverse  of  cumulative  distribution  function 

Importance  sampling 

Load  or  length 

Load  and  resistance  factor  design 

Monte  Carlo  simulation 

Moments  (mean  and  standard  deviation) 

Number  of  simulation  cycles 

Probability 

Parameters 

Probability  density  function 
Strength  or  resistance 

Software  for  reliability  and  stability  assessment  of  concrete 
gravity  structures 
Safety  factor 
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C1 


Stddv 

t 

Var 

X 

Y 
Z 
a 

P 

A 

e 

Y 
8 
4> 

<t) 

X 

\i 

o 

P 

a 


Standard  deviation 

Variance 
Random  variable 
Random  variable 
Performance  function 
Directional  cosine 
Reliability  index 
Very  small  quantity 
Standard  error  for  regression 
Partial  safety  factor 
Tolerance 

Cumulative  distribution  function  of  standard  normal  variate 
Probability  density  function  of  standard  normal  variate 
Eigenvalue 
Mean 

Standard  deviation 
Correlation  coefficient 
Partial  derivative 


Subscripts 

i  f'  iteration  (or  component) 

m  Mean 

s  Shifted 

t  Target 

u  Unsatisfactory  performance 

X  Random  variable 


Superscript 

N  Equivalent  normal  distribution 

*  Design  point 

-1  Inverse  of  a  function 
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Analysis  of  Plane  Frame  Structures  (CFRAME) 

Jan  1983 

Instruction  Report  K-83*‘2 

User’s  Guide:  Computer  Program  for  Generation  of  Engineering 
Geometry  (SKETCH) 

Jun  1983 

Instruction  Report  K-83'5 

User’s  Guide:  Computer  Program  to  Calculate  Shear,  Moment, 
and  Thrust  (CSMT)  from  Stress  Results  of  a  Two-Dimensional 

Finite  Element  Analysis 

Jul  1983 

Technical  Report  K-83-1 

Basic  Pile  Group  Behavior 

Sep  1983 

Technical  Report  K-83-3 

Reference  Manual:  Computer  Graphics  Program  for  Generation  of 
Engineering  Geometry  (SKETCH) 

Sep  1983 

Technical  Report  K-83-4 

Case  Study  of  Six  Major  General-Purpose  Finite  Element  Programs 

Oct  1983 

Instruction  Report  K-84-2 

User’s  Guide:  Computer  Program  for  Optimum  Dynamic  Design 
of  Nonlinear  Metal  Plates  Under  Blast  Loading  (CSDOOR) 

Jan  1984 

Instruction  Report  K-84-7 

User’s  Guide:  Computer  Program  for  Determining  Induced 

Stresses  and  Consolidation  Settlements  (CSETT) 

Aug  1984 

Instrucrtlon  Report  K-84-8 

Seepage  Analysis  of  Confined  Flow  Problems  by  the  Method  of 
Fragments  (CFRAG) 

Sep  1984 

Instruction  Report  K-84-1 1 

User’s  Guide  for  Computer  Program  CGFAG,  Concrete  General 

Flexure  Analysis  with  Graphics 

Sep  1984 

Technical  Report  K-84-3 

Computer-Aided  Drafting  and  Design  for  Corps  Structural 

Engineers 

Oct  1984 

Technical  Report  ATC-86-5 

Decision  Logic  Table  Formulation  of  ACI  318-77,  Building  Code 
Requirements  for  Reinforced  Concrete  for  Automated  Con¬ 
straint  Processing,  Volumes  1  and  II 

Jun  1986 

Technical  Report  ITL-87-2 

A  Case  Committee  Study  of  Finite  Element  Analysis  of  Concrete 

Flat  Slabs 

Jan  1987 

Instruction  Report  lTL-87-1 

User’s  Guide:  Computer  Program  for  Two-Dimensional  Analysis 
of  U-Frame  Structures  (CUFRAM) 

Apr  1987 

Instruction  Report  ITL-87-2 

User’s  Guide:  For  Concrete  Strength  Investigation  and  Design 
(CASTR)  in  Accordance  with  ACI  318-83 

May  1987 

Technical  Report  ITL-87-6 

Finite-Element  Method  Package  for  Solving  Steady-State  Seepage 
Problems 

May  1987 

Instruction  Report  ITL-87-3 

User’s  Guide:  A  Three-Dimensional  Stability  Analysis/Design 

Program  (3DSAD)  Module 

Report  1 :  Revision  1 :  General  Geometry 

Report  2:  General  Loads  Module 

Report  6:  Free-Body  Module 

Jun  1987 

Jun  1987 
Sep  1989 
Sep  1989 
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Instruction  Report  ITL-87-4 

User’s  Gui(je:  2-D  Frame  Analysis  Link  Program  (L1NK2D) 

Jun  1987 

Technical  Report  lTL-87-4 

Finite  Element  Stuidies  of  a  Horizontally  Framed  Miter  Gate 

Report  1 :  Initial  and  Refined  Finite  Element  Models  (Phases 

A,  B,  and  C),  Volumes  1  and  II 

Report  2:  Simplified  Frame  Model  (Phase  D) 

Report  3:  Alternate  Configuration  Miter  Gate  Finite  Element 
Studies-Open  Section 

Report  4:  Alternate  Configuration  Miter  Gate  Finite  Element 
Studies-Closed  Sections 

Report  5:  Alternate  Configuration  Miter  Gate  Finite  Element 
Studies-Additional  Closed  Sections 

Report  6:  Elastic  Buckling  of  Girders  In  Horizontally  Framed 

Miter  Gates 

Report  7:  Application  and  Summary 

Aug  1987 

Instruction  Report  GL-87-1 

User’s  Guide:  UTEXAS2  Slope-Stability  Package;  Volume  1, 

User’s  Manual 

Aug  1987 

Instruction  Report  ITL-87-5 

Sliding  Stability  of  Concrete  Structures  (CSLIDE) 

Oct  1987 

Instruction  Report  ITL-87-6 

Criteria  Specifications  for  and  Validation  of  a  Computer  Program 
for  the  Design  or  Investigation  of  Horizontally  Framed  Miter 

Gates  (CMITER) 

Dec  1987 

Technical  Report  ITL-87-8 

Procedure  for  Static  Analysis  of  Gravity  Dams  Using  the  Finite 

Element  Method  -  Phase  1  a 

Jan  1 988 

Instruction  Report  ITL-88-1 

User’s  Guide:  Computer  Program  for  Analysis  of  Planar  Grid 
Structures  (CGRID) 

Feb  1988 

Technical  Report  ITL-88-1 

Development  of  Design  Formulas  for  Ribbed  Mat  Foundations 
on  Expansive  Soils 

Apr  1988 

Technical  Report  ITL-88-2 

User’s  Guide:  Pile  Group  Graphics  Display  (CPGG)  Post¬ 
processor  to  CPGA  Program 

Apr  1988 

Instruction  Report  ITL-88-2 

User’s  Guide  for  Design  and  Investigation  of  Horizontally  Framed 

Miter  Gates  (CMITER) 

Jun  1988 

Instruction  Report  ITL-88-4 

User’s  Guide  for  Revised  Computer  Program  to  Calculate  Shear, 
Moment,  and  Thrust  (CSMT) 

Sep  1 988 

Instruction  Report  GL-87-1 

User’s  Guide:  UTEXAS2  Slope-Stability  Package;  Volume  II, 

Theory 

Feb  1989 

Technical  Report  ITL-89-3 

User’s  Guide:  Pile  Group  Analysis  (CPGA)  Computer  Group 

Jul  1989 

Technical  Report  ITL-89-4 

CBASIN-Structural  Design  of  Saint  Anthony  Falls  Stilling  Basins 
According  to  Corps  of  Engineers  Criteria  for  Hydraulic 

Structures;  Computer  Program  X0098 

Aug  1 989 
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Technical  Report  ITL-89-5 

CCHAN-Structural  Design  of  Rectangular  Channels  Accor(jing 
to  Corps  of  Engineers  Criteria  for  Hydraulic 

Structures;  Computer  Program  X0097 

Aug  1 989 

Technical  Report  ITL*89-6 

The  Response-Specitrum  Dynamic  Analysis  of  Gravity  Dams  Using 
the  Finite  Element  Method;  Phase  II 

Aug  1 989 

Contract  Report  ITL-89-1 

State  of  the  Art  on  Expert  Systems  Applications  in  Design, 

Construction,  and  Maintenance  of  Structures 

Sep  1 989 

Instruction  Report  lTL-90-1 

User’s  Guide;  Computer  Program  for  Design  and  Analysis 
of  Sheet  Pile  Walls  by  Classical  Methods  (CWALSHT) 

Feb  1990 

Technical  Report  ITL-90-3 

Investigation  and  Design  of  U-Frame  Structures  Using 

Program  CUFRBC 

Volume  A:  Program  Criteria  and  Documentation 

Volume  B:  User’s  Guide  for  Basins 

Volume  C:  User’s  Guide  for  Channels 

May  1990 

Instruction  Report  ITL-90-6 

User’s  Guide:  Computer  Program  for  Two-Dimensional  Analysis 
of  U-Frame  or  W- Frame  Structures  (CWFRAM) 

Sep  1 990 

Instruction  Report  ITL-90-2 

User’s  Guide:  Pile  Group-Concrete  Pile  Analysis  Program 
(CPGC)  Preprocessor  to  CPGA  Program 

Jun  1990 

Technical  Report  ITL-91-3 

Application  of  Finite  Element,  Grid  Generation,  and  Scientific 
Visualization  Techniques  to  2-D  and  3-D  Seepage  and 

Groundwater  Modeling 

Sep  1990 

Instruction  Report  ITL-91-1 

User’s  Guide:  Computer  Program  for  Design  and  Analysis 
of  Sheet-Pile  Walls  by  Classical  Methods  (CWALSHT) 

Including  Rowe’s  Moment  Reduction 

Oct  1991 

Instruction  Report  ITL-87-2 
(Revise(d) 

User’s  Guide  for  Concrete  Strength  Investigation  and  Design 
(CASTR)  in  Accordance  with  ACI  318-89 

Mar  1992 

Technical  Report  ITL-92-2 

Finite  Element  Modeling  of  Welded  Thick  Plates  for  Bonneville 
Navigation  Lock 

May  1992 

Technical  Report  lTL-92-4 

Introduction  to  the  Computation  of  Response  Spectrum  for 

Earthquake  Loading 

Jun  1992 

Instruction  Report  ITL-92-3 

Concept  Design  Example,  Computer-Aided  Structural 

Modeling  (CASM) 

Report  1 :  Scheme  A 

Report  2:  Scheme  B 

Reports:  Scheme  C 

Jun  1992 
Jun  1992 
Jun  1992 

Instruction  Report  ITL-92-4 

User’s  Guide:  Computer-Aided  Structural  Modeling 
(CASM)  -Version  3.00 

Apr  1992 

Instruction  Report  ITL-92-5 

Tutorial  Guide:  Computer-Aided  Structural  Modeling 
(CASM)  -Version  3.00 

Apr  1992 
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Contract  Report  ITL-92-1 

Optimization  of  Steel  Pile  Foundations  Using  Optimality  Criteria 

Jun  1992 

Technical  Report  ITL~92-7 

Refined  Stress  Analysis  of  Melvin  Price  Locks  and  Dam 

Sep  1992 

Contract  Report  ITL-92-2 

Knowledge-Based  Expert  System  for  Selection  and  Design 
of  Retaining  Structures 

Sep  1992 

Contract  Report  ITL-92-3 

Evaluation  of  Thermal  and  Incremental  Construc^tion  Effects 
for  Monoliths  AL-3  and  AL-5  of  the  Melvin  Price  Locks 
and  Dam 

Sep  1992 

Instruction  Report  GL-87-1 

User’s  Guide:  UTEXAS3  Slope-Stability  Package;  Volume  IV, 

User’s  Manual 

Nov  1992 

Technical  Report  ITL-92-1 1 

The  Seismic  Design  of  Waterfront  Retaining  Structures 

Nov  1992 

Technical  Report  ITL-92-1 2 

Computer-Aided,  Field-Verified  Structural  Evaluation 

Report  1 :  Development  of  Computer  Modeling  Techniques 
for  Miter  Lock  Gates 

Nov  1992 

Report  2:  Field  Test  and  Analysis  Correlation  at  John  Hollis 
Bankhead  Lock  and  Dam 

Dec  1992 

Report  3:  Field  Test  and  Analysis  Correlation  of  a  Vertically 

Framed  Miter  Gate  at  Emsworth  Lock  and  Dam 

Dec  1993 

Instruction  Report  GL-87-1 

User’s  Guide:  UTEXAS3  Slope-Stability  Package;  Volume  III, 

Example  Problems 

Dec  1992 

Technical  Report  ITL-93-1 

Theoretical  Manual  for  Analysis  of  Arch  Dams 

Jul  1993 

Technical  Report  ITL-93-2 

Steel  Structures  for  Civil  Works,  General  Considerations 
for  Design  and  Rehabilitation 

Aug  1993 

Technical  Report  ITL-93-3 

Soil-Structure  Interaction  Study  of  Red  River  Lock  and  Dam 

No.  1  Subjected  to  Sediment  Loading 

Sep  1993 

Instruction  Report  ITL-93-3 

User’s  Manual— ADAP,  Graphics-Based  Dam  Analysis  Program 

Aug  1993 

Instruction  Report  ITL-93-4 

Load  and  Resistance  Factor  Design  for  Steel  Miter  Gates 

Oct  1993 

Technical  Report  ITL-94-2 

User’s  Guide  for  the  Incremental  Construction,  Soil-Structure  Interaction 
Program  SOILSTRUCT  with  Far-Field  Boundary  Elements 

Mar  1994 

Instruction  Report  ITL-94-1 

Tutorial  Guide:  Computer-Aided  Structural  Modeling  (CASM); 

Version  5.00 

Apr  1994 

Instruction  Report  ITL-94-2 

User’s  Guide:  Computer-Aided  Structural  Modeling  (CASM); 

Version  5.00 

Apr  1994 

Technical  Report  ITL-94-4 

Dynamics  of  Intake  Towers  and  Other  MDOF  Structures  Under 
Earthquake  Loads:  A  Computer-Aided  Approach 

Jul  1994 

Technical  Report  ITL-94-5 

Procedure  for  Static  Analysis  of  Gravity  Dams  Including  Foundation 
Effects  Using  the  Finite  Element  Method  -  Phase  1 B 

Jul  1994 
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Instruction  Report  ITL-94-5 

User’s  Guide:  Computer  Program  for  Winkler  Soil-Structure 

Interaction  Analysis  of  Sheet-Pile  Walls  (CWALSSI) 

Nov  1994 

Instruction  Report  ITL>94-6 

User’s  Guide:  Computer  Program  for  Analysis  of  Beam-Column 
Structures  with  Nonlinear  Supports  (CBEAMC) 

Nov  1994 

Instruction  Report  ITL-94-7 

User’s  Guide  to  CTWALL  -  A  Microcomputer  Program  for  the 

Analysis  of  Retaining  and  Flood  Walls 

Dec  1994 

Contract  Report  ITL-95-1 

Comparison  of  Barge  Impact  Experimental  and  Finite  Element 

Results  for  the  Lower  Miter  Gate  of  Lock  and  Dam  26 

Jun  1995 

Technical  Report  ITL-95-5 

Soil-Structure  Interaction  Parameters  for  Structured/Cemented 

Silts 

Aug  1995 

instruction  Report  ITL~95-1 

User’s  Guide:  Computer  Program  for  the  Design  and  Investigation 
of  Horizontally  Framed  Miter  Gates  Using  the  Load  and  Resistance 
Factor  Criteria  (CMITER-LRFD) 

Aug  1995 

Technical  Report  ITL-95'8 

Constitutive  Modeling  of  Concrete  for  Massive  Concrete  Structures, 

A  Simplified  Overview 

Sep  1995 

Instruction  Report  ITL-96-1 

User’s  Guide:  Computer  Program  for  Two-Dimensional  Dynamic 
Analysis  of  U-Frame  or  W-Frame  Structures  (CDWFRM) 

Jun  1996 

Instruction  Report  ITL-96-2 

Computer-Aided  Structural  Modeling  (CASM),  Version  6.00 

Report  1 :  Tutorial  Guide 

Report  2:  User’s  Guide 

Report  3:  Scheme  A 

Report  4:  Scheme  B 

Report  5:  Scheme  C 

Jun  1996 

Technical  Report  ITL-96-8 

Hyperbolic  Stress-Strain  Parameters  for  Structured/Cemented  Silts 

Aug  1996 

Instruction  Report  ITL-96-3 

User’s  Guide:  Computer  Program  for  the  Design  and  Investigation 
of  Horizontally  Framed  Miter  Gates  Using  the  Load  and  Resistance 
Factor  Criteria  (CMITERW-LRFD)  Windows  Version 

Sep  1996 

Instruction  Report  ITL-97-1 

User’s  Guide:  Computer  Aided  Inspection  Forms  for  Hydraulic  Steel 
Structures  (CAIF-HSS),  Windows  Version 

Sep  1997 

Instrucrtion  Report  ITL-97-2 

User’s  Guide:  Arch  Dam  Stress  Analysis  System  (ADSAS) 

Aug  1997 

Instruction  Report  ITL-98-1 

User’s  Guide  for  the  Three-Dimensional  Stability  Analysis/Design 
(3DSAD)  Program 

Sep  1998 

Technical  Report  ITL-98-4 

Investigation  of  At-Rest  Soil  Pressures  due  to  Irregular  Sloping  Soil 
Surfaces  and  CSOILP  User’s  Guide 

Sep  1998 

Technical  Report  ITL-98'5 

The  Shear  Ring  Method  and  the  Program  Ring  Wall 

Sep  1998 

Technical  Report  ITL-98‘6 

Reliability  and  Stability  Assessment  of  Concrete  Gravity  Structures 
(RCSLIDE):  Theoretical  Manual 

Dec  1998 

